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[i]  Despite  evidence  from  ground-based  data  that  flow  over  mountains  is  a  dominant 
source  of  gravity  waves  (GWs)  for  the  Northern  Hemisphere  winter  middle  atmosphere, 

GW-related  signals  in  global  limb  radiances  from  the  Microwave  Limb  Sounder  (MLS)  on 
the  Upper  Atmosphere  Research  Satellite  (UARS)  have  shown  little  direct  evidence  of 
mountain  waves.  We  address  this  issue  by  combining  a  renewed  analysis  of  MLS  limb-track 
and  limb-scan  radiances  with  global  mountain  wave  modeling  using  the  Naval  Research 
Laboratory  Mountain  Wave  Forecast  Model  (MWFM).  MLS  radiance  variances  show 
characteristics  consistent  with  mountain  waves,  such  as  enhanced  variance  over  specific 
mountain  ranges  and  annual  variations  that  peak  strongly  in  winter.  However,  direct 
comparisons  of  MLS  variance  maps  with  MWFM-simulated  mountain  wave  climatologies 
reveal  limited  agreement.  We  further  develop  a  detailed  “MLS  GW  visibility  function”  that 
accurately  specifies  the  three-dimensional  in-orbit  sensitivity  of  the  MLS  limb-track 
radiance  measurement  to  a  spectrum  of  GWs  with  different  wavelengths  and  horizontal 
propagation  directions.  On  postprocessing  MWFM-generated  mountain  wave  fields 
through  these  MLS  visibility  filters,  we  generate  MWFM  variance  maps  that  agree 
substantially  better  with  MLS  radiance  variances.  This  combined  data  analysis  and  MLS- 
filtered  MWFM  modeling  leads  us  to  conclude  that  many  MLS  variance  enhancements  can 
be  associated  with  mountain  waves  forced  by  flow  over  specific  mountainous  terrain.  These 
include  mountain  ranges  in  Europe  (e.g.,  Scandinavia;  Alps;  Scotland;  Ural,  Putoran, 

Altai,  Hangay  and  Sayan  Mountains;  Yablonovyy,  Stanovoy,  Khingan,  Verkhoyansk  and 
Central  Ranges),  North  America  (e.g.,  Brooks  Range,  MacKenzie  Mountains,  Colorado 
Rockies),  southeastern  Greenland,  and  Iceland.  Our  results  show  that  given  careful 
consideration  of  the  in-orbit  sensitivity  of  the  instrument  to  GWs,  middle  atmospheric  limb 
radiances  measured  from  UARS  MLS,  as  well  as  from  the  new  MLS  instrument  on  the  Earth 
Observing  System  (EOS)  satellite,  can  provide  important  global  information  on  mountain 
waves  in  the  extratropical  Northern  Hemisphere  stratosphere  and  mesosphere.  Index 
TERMS:  3334  Meteorology  and  Atmospheric  Dynamics:  Middle  atmosphere  dynamics  (0341,  0342);  3360 
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1.  Introduction 

[2]  Atmospheric  gravity  waves  (GWs)  play  cmcial  roles 
in  driving  the  general  circulation  and  thermal  structure  of  the 
atmosphere.  Quasi-continuous  GW  breaking  at  a  range  of 
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heights  around  the  globe  maintains  body  forces  and  turbu¬ 
lent  diffusion  that  drive  winds,  temperatures  and  chemical 
constituent  distributions  well  away  from  those  of  an  atmo¬ 
sphere  lacking  these  processes  [e.g.,  Andrews ,  1987;  Holton 
and  Alexander ,  2000;  McIntyre ,  2001;  Fritts  and  Alexander , 
2003].  Since  the  spectrum  of  breaking  gravity  waves  is  not 
resolved  in  global  climate  and  weather  prediction  models  at 
present,  these  unresolved  GW  processes  must  be  compre¬ 
hensively  parameterized  so  that  models  can  reproduce  real¬ 
istic  circulations  on  both  short  and  long  timescales  [see,  e.g., 
Hamilton ,  1996;  McLandress ,  1998;  Kim  et  al. ,  2003]. 
Gravity  waves  can  produce  other  important  effects,  such  as 
wave  clouds  at  various  heights  that  can  have  important 
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follow-on  effects  on  processes  such  as  precipitation  and 
ozone  [e.g.,  Jensen  and  Toon ,  1994;  Carslaw  et  al,  1998; 
Koch  and  Siedlarz ,  1999;  Thayer  et  al ,  2003].  This  has  led 
to  a  new  class  of  parameterizations  of  subgrid-scale  gravity 
wave  effects  on  clouds  for  global  models  [e.g.,  Cusack  et  al , 
1999;  Bacmeister  et  al,  1999;  Pierce  et  al,  2003]. 

[3]  GWs  are  generated  by  a  variety  of  processes  in  the 
lower  atmosphere  [ Fritts  and  Alexander ,  2003].  Of  these, 
flow  over  mountains  is  believed  to  be  one  of  the  dominant 
sources,  particularly  in  the  extratropics  during  winter  [e.g., 
Nastrom  and  Fritts ,  1992].  The  drag  produced  by  the 
breaking  of  these  so-called  mountain  waves  is  important 
for  both  short-term  weather  prediction  and  long-term  cli¬ 
mate  modeling,  and  so  considerable  effort  has  been  devoted 
to  mountain  wave  drag  parameterization  schemes  for  global 
models  [see,  e.g.,  Kim  et  al,  2003]. 

[4]  Unfortunately,  many  aspects  of  GW  parameterizations 
are  uncertain  and  highly  simplified,  mainly  because  of  a 
lack  of  detailed  global  data  on  GWs  that  could  be  used  to 
intercompare,  assess  and  constrain  the  various  models  of 
these  processes  currently  in  use.  This  paucity  of  global  data 
has  arisen  because,  until  recently,  limb  and  nadir  viewing 
remote-sensing  instruments  on  satellites  lacked  the  spatial 
or  temporal  resolution  to  resolve  atmospheric  GW  pertur¬ 
bations.  Fetzer  and  Gille  [1994],  however,  showed  that  raw 
temperature  data  acquired  by  the  Limb  Infrared  Monitor  of 
the  Stratosphere  (LIMS)  on  the  Nimbus-7  satellite  resolved 
long  wavelength  gravity  wave  fluctuations.  During  the  last 
decade,  other  high-resolution  satellite  instruments  have 
attained  the  necessary  resolutions  and  accuracies  to  provide 
our  first  tentative  global  glimpses  of  these  long  wavelength 
portions  of  the  GW  spectrum.  One  of  the  most  important 
instruments  in  this  critical  emerging  database  has  been  the 
Microwave  Limb  Sounder  (MLS)  on  the  Upper  Atmosphere 
Research  Satellite  (UARS)  [ Waters  et  al,  1999].  Since  the 
mid-1990s,  UARS  MLS  63  GHz  radiance  fluctuations  have 
been  used  to  study  global  GW  morphologies  in  the  strato¬ 
sphere  and  mesosphere  [Wu  and  Waters ,  1996a,  1996b, 
1997;  Alexander ,  1998;  McLandress  et  al,  2000;  Wu ,  2001; 
Jiang  and  Wu ,  2001;  Wu  and  Jiang ,  2002;  Jiang  et  al, 
2002,  2003,  2004]. 

[5]  Previous  global  maps  of  MLS  radiance  variances  in 
the  stratosphere  and  mesosphere  showed  good  correlations 
between  enhanced  variance  and  strong  background 
wind  speeds  [e.g.,  Wu  and  Waters ,  1996a,  1996b,  1997; 
McLandress  et  al,  2000;  Jiang  and  Wu ,  2001].  Alexander 
[1998]  showed  that  much  of  this  structure  could  be 
reproduced  by  a  simple  source-independent  GW  ray¬ 
tracing  model,  which  accounted  for  the  “Doppler-shifting” 
effects  of  background  winds  on  the  gravity  wave  spectrum 
and  the  finite  vertical  width  of  the  MLS  channel  weighting 
functions,  which  allowed  only  those  waves  with  long 
vertical  wavelengths  to  be  resolved.  McLandress  et  al 
[2000]  performed  a  detailed  follow-up  study  of  MLS 
variances  measured  by  channels  3  and  13  (altitudes  z  ~ 
38  km)  and  used  a  similar  global  GW  ray  model  to  simulate 
the  observed  variance  distributions.  Their  work  showed  that 
refraction  and  observational  filtering  of  waves  by  back¬ 
ground  winds  alone  could  not  explain  all  of  the  structure  in 
MLS  GW  variance  maps.  For  example,  model  sensitivity 
experiments  suggested  that  large  increases  in  stratospheric 
variances  in  the  summer  subtropics  and  in  winter  over 


the  southern  Andes  were  consistent  with  intense  local 
generation  of  gravity  waves  from  underlying  convective 
and  topographic  sources,  respectively. 

[6]  Eckermann  and  Preusse  [1999]  analyzed  a  week’s 
worth  of  infrared  limb  temperature  data  acquired  by  the 
CRISTA  instrument  on  the  Shuttle  Pallet  Satellite  in 
November  1994.  After  isolating  stratospheric  temperature 
perturbations,  they  used  simple  theory  and  output  from  the 
Naval  Research  Laboratory  (NRL)  Mountain  Wave  Forecast 
Model  (MWFM)  to  show  that  enhanced  stratospheric  tem¬ 
perature  fluctuations  over  the  southern  tip  of  South  America 
were  produced  by  mountain  waves  forced  by  flow  across  the 
Andes.  A  detailed  follow-up  study  by  Preusse  et  al  [2002] 
confirmed  this  using  retrieval  modeling,  GW  theory,  a 
mesoscale  model  and  ray-tracing  experiments.  Given  these 
indications,  Jiang  et  al  [2002]  conducted  a  thorough  study  of 
stratospheric  radiance  variances  over  the  southern  Andes 
region,  utilizing  slightly  more  than  2  years  of  MLS  limb 
track  data  between  late  1994  and  early  1997.  Variances  over 
the  Andes  exhibited  a  strong  annual  variation,  peaking  in 
winter,  with  somewhat  weaker  wintertime  activity  in  1996 
compared  to  1995.  These  features  were  reproduced  quite  well 
by  a  detailed  MWFM  “hindcasf  ’  simulation  formulated  for 
this  particular  time  period.  Wu  and  Jiang  [2002]  also  reported 
enhancements  of  MLS  variances  over  Antarctica  that  may 
also  be  mountain  wave-related. 

[7]  Mountain  waves  are  believed  to  be  more  energetic  and 
geographically  prevalent  in  the  Northern  Hemisphere 
[. Bacmeister ,  1993;  Fritts  and  Alexander ,  2003],  because  of 
the  more  extended  mountainous  landmasses  in  this  hemi¬ 
sphere.  Yet  little  if  any  progress  has  been  made  to  date  in 
extracting  unambiguous  data  on  mountain  waves  from  MLS 
radiances  in  the  Northern  Hemisphere.  McLandress  et  al 
[2000]  documented  highly  structured  variance  distributions 
in  the  northern  extratropics,  but  they  were  unable  to  draw 
conclusions  as  to  the  sources  of  this  wave  activity  because 
strong  longitudinal  and  day-to-day  variations  in  stratospheric 
wind  fields  complicated  their  model-data  comparisons.  Jiang 
and  Wu  [2001]  further  demonstrated  how  strong  variable 
stratospheric  flow  fields  in  the  north  substantially  modulated 
MLS  variances,  seemingly  masking  any  signatures  of  wave 
sources.  Such  difficulties  do  not  appear  to  be  unique  to  MLS. 
Despite  showing  clear  evidence  of  enhanced  stratospheric 
temperature  variance  from  gravity  waves  over  deep  tropical 
convection,  GP S/MET  occultation  data  have  to  date  yielded 
little  clear  evidence  of  enhanced  temperature  variances  over 
mountains  in  the  Northern  Hemisphere  during  winter  [Tsuda 
et  al,  2000]. 

[8]  Despite  these  complications,  MLS  satellite  data  do 
show  occasional  hints  of  variance  enhancements  over  topog¬ 
raphy  during  northern  winter.  For  example,  Figure  5b  of 
McLandress  et  al  [2000]  shows  enhancements  in  limb  track 
variances  over  mountainous  regions  such  as  Scandinavia, 
central  Eurasia,  Ellesmere  Island  and  southern  Greenland, 
where  long  wavelength  stratospheric  mountain  waves  have 
sometimes  been  observed  [e.g.,  Leutbecher  and  Volkert , 
2000;  Dornbrack  et  al,  2002].  Furthermore,  Eckermann 
and  Preusse  [1999]  found  enhanced  CRISTA  temperature 
perturbations  over  central  Eurasia  and  used  MWFM  hind- 
casts  to  show  that  they  were  consistent  with  energetic 
stratospheric  mountain  waves  emanating  from  flow  across 
significant  underlying  mountain  ranges  in  this  region. 
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[9]  Why  are  mountain  waves  evident  in  some  satellite 
observations,  yet  seemingly  absent  or  obscured  in  others? 
Do  MLS  radiances  contain  any  information  on  stratospheric 
mountain  waves  in  the  Northern  Hemisphere?  If  not,  why 
not,  particularly  when  clear  mountain  wave  signals  appear 
in  MLS  data  from  the  Southern  Hemisphere  [Jiang  et  al ., 
2002]?  If  there  is  mountain  wave  information  buried  within 
these  data,  can  we  extract  it  to  provide  much-needed  global 
information  on  these  waves?  Motivated  by  these  questions, 
this  paper  combines  a  careful  analysis  of  MLS  stratospheric 
radiance  data  from  the  Northern  Hemisphere  with  detailed 
companion  modeling  of  anticipated  global  distributions  of 
mountain  waves  that  MLS  might  resolve.  Our  scientific 
approach  is  as  follows:  (1)  construct  improved  climatolog¬ 
ical  maps  of  GW-related  MLS  variances  in  the  Northern 
Hemisphere;  (2)  identify  specific  zones  of  variance  en¬ 
hancement  that  are  the  likeliest  candidates  for  a  mountain 
wave  explanation;  (3)  analyze  variances  over  these  “focus 
regions”  using  all  available  MLS  data  in  ways  that  test  the 
mountain  wave  hypothesis  further;  (4)  conduct  detailed 
“hindcasf  ’  modeling  of  these  observations  using  the  NRL 
MWFM,  paying  particular  attention  to  the  instrumental  and 
environmental  effects  that  influence  how  mountain  waves 
manifest  in  MLS  radiances. 

[10]  In  section  2,  we  describe  the  MLS  GW  variance 
calculation  techniques  and  their  recent  improvements.  In 
section  3,  we  analyze  MLS  radiance  variances  in  Northern 
Hemispheric  winter  and  identify  geographical  regions  with 
potential  mountain  wave-induced  variance  enhancements. 
In  section  4,  we  outline  our  modeling  strategy,  which  uses 
the  MWFM  to  simulate  MWs  that  manifest  in  MLS 
radiances.  In  particular,  we  develop  a  complete  three- 
dimensional  (3-D)  MLS  visibility  function  that  mimics 
the  MLS  GW  observation,  specifically  its  sensitivity  to 
waves  of  different  wavelengths  and  horizontal  orientations 
with  respect  to  the  MLS  line-of-slight  (LOS).  Section  5 
compares  modeled  MWFM  fields  with  the  MLS  variance 
data.  Results  are  briefly  summarized  and  discussed  in 
section  6. 

2.  Computation  of  MLS  Radiance  Variances 

[11]  The  UARS  MLS  consists  of  three  double-sideband 
radiometers  that  measure  atmospheric  02,  03,  CIO,  and 
H20  emission  features  near  63,  183  and  205  GHz  [ Waters , 
1993;  Barath  et  al ,  1993].  The  63 -GHz  radiometer  has 
15  spectral  channels,  with  channel  8  measuring  the  central 
02  line,  and  progressing  to  channels  1/15  that  are  located  at 
±200  MHz  from  this  central  line.  In  limb-scan  operation 
mode  (years  1991-1997),  MLS  step-scans  the  atmospheric 
limb  between  90  km  and  the  surface  in  65  s  with  ^2-s 
integration  time  for  each  measurement.  The  limb  radiances 
become  saturated  as  the  antenna  views  tangent  heights  near 
the  surface.  In  limb-track  operation  mode  (years  1994- 
1997),  the  MLS  antenna  tracks  a  constant  tangent  height 
(often  at  ^18  km)  where  the  limb  radiances  are  saturated. 
UARS  moves  at  a  constant  speed  of  7.5  km  s-1  along  its 
orbit  and  the  MLS  antenna  points  orthogonal  to  the  satellite 
velocity.  Thus  the  horizontal  separation  between  adjacent 
limb-track  measurements  is  ~15  km. 

[12]  The  MLS  saturated  radiance  basically  measures  the 
atmospheric  temperature  of  the  saturation  layer,  and  wave- 


induced  temperature  fluctuations  can  appear  as  along-track 
fluctuations  in  measured  radiances  [Wu  and  Waters ,  1996b]. 
To  derive  GW  variances  from  MLS  saturated  radiances,  Wu 
and  Waters  [1996b]  used  a  straightforward  analysis  method 
consisting  of  three  main  parts.  Our  current  version  of  that 
basic  analysis  method  is  outlined  below. 

2.1.  Variance  Estimate 

[13]  For  each  MLS  63GHz  frequency  channel,  the  esti¬ 
mated  radiance  variance  b2  is  obtained  from  n  consecutive 
individual  measurements  of  saturated  radiances,  as  follows: 

52  =  ST^i-a-bzj)2,  (la) 

i=  1 

where  yt  and  zt  are  individual  radiance  and  tangent  height 
measurements,  respectively.  Wu  and  Waters  [1996a,  1996b] 
used  n  =  6  for  both  limb-scan  and  limb-track  variances.  In 
this  study,  we  use  n  =  4  for  limb-scan  and  n  =  6  for  limb- 
track  data,  for  reasons  described  in  section  2.4.  Because  the 
saturated  radiances  may  depend  slightly  on  tangent  height 
as  the  limb  path  length  decreases,  a  linear  trend  (described 
by  parameters  a  and  b )  is  fitted  to  the  radiances  to  remove 
both  the  mean  radiance  trend  and  any  tangent  height 
dependence.  The  factor  n  —  2  in  the  denominator  of 
equation  (la)  comes  from  reduced  degrees  of  freedom  in 
this  variance  estimate  [Wu  and  Waters ,  1997],  since  a  and  b 
are  two  constraint  parameters  for  yt  and  zt. 

2.2.  Variance  Averaging 

[14]  The  next  step  in  the  analysis  involves  substantial 
averaging  of  many  individual  variance  estimates  aj 

d2=p-'J2°j  (ib) 

j= 1 

[15]  Since  the  small-scale  atmospheric  (wave-induced) 
radiance  variances  are  often  transient  and  weak  compared 
to  instrument  noise  [e.g.,  McLandress  et  al ,  2000;  Wu , 
2001],  the  statistical  uncertainties  in  the  estimated  varian¬ 
ces  must  be  reduced  via  substantial  averaging  as  in 
equation  (lb)  (p  1).  To  a  good  approximation,  the 
statistical  uncertainty  (random  error)  of  one  given  MLS 
radiance  variance  estimate  a2  can  be  expressed  as  [Wu 
and  Waters ,  1997] 

|  o’2  —  cj2  |  «  v^cr2,  (2) 

where  cr2  is  the  true  variance.  This  measurement  uncertainty 
on  the  right-hand  side  of  equation  (2)  is  reduced  by  a  factor 
of  y/p  when  /^-independent  variance  estimates  are  averaged 
as  in  equation  (lb)  [Wu  and  Waters ,  1997].  Such  averaging 
is  achieved  at  the  expense  of  the  final  spatial  and  temporal 
resolution  of  the  variance  maps,  and  thus  some  tradeoff 
between  signal  averaging  and  final  resolution  must  be 
made.  In  practice,  improvements  from  the  averaging  are 
also  ultimately  limited  by  other  factors  such  as  the 
atmospheric/noise  variance  ratio,  spatial/temporal  resolu¬ 
tion,  and  MLS  sampling  rate.  The  statistical  uncertainty  in 
the  mean  variance,  discussed  here,  should  not  be  confused 
with  the  large  standard  deviations  of  the  MLS  variances  that 
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Table  1.  Estimated  MLS  63  GHz  Channel  Noise  Variances,  for 
Limb-Track  (Third  Column)  and  Limb-Scan  (Fourth  Column) 
Variances,  Compared  With  the  Calibrated  Noise  Values  From  Wu 
and  Waters  [1991  f 


Combined 

Channels 

Approximate 
Height,  km 

Limb-Track 

#(K2) 

Limb -Scan 

i  (k2) 

(Noise)2 

2 

1/15 

28 

0.0018 

0.0022 

0.0025 

2/14 

33 

0.0027 

0.0034 

0.0032 

3/13 

38 

0.0048 

0.0051 

0.0072 

4/12 

43 

0.0096 

0.010 

0.016 

5 

48 

0.038 

0.039 

0.034 

6/10 

53 

0.039 

0.040 

0.068 

7/9 

61 

0.070 

0.071 

0.12 

8 

80 

0.13 

0.13 

0.10 

aNote  that  most  channel  data  are  combined,  except  channels  5  and  11,  to 
reduce  the  net  noise  variance  (channel  1 1  did  not  function  during  the  UARS 
MLS  mission).  The  noise  values  in  the  last  column  are  from  Wu  and  Waters 
[1997]. 

can  arise  because  of  geophysical  variations  in  the  intensity 
of  the  local  wave  fields  (see,  e.g.,  Figure  7  of  McLandress  et 
al.  [2000]). 

2.3.  Noise  Estimation  and  Removal 

[i6]  The  final  step  in  the  process  involves  estimating 
and  removing  instrument  noise  variances.  The  variance  in 
equation  (la)  has  two  main  components:  instrument  noise 
and  atmospheric  variance,  i.e.,  a2  =  cr^  +  ga  (here  we 
neglect  an  additional  nonlinear  pressure  error  variance 
c i2nl  that  arises  only  for  data  from  channels  1/15  and 
2/14  [Wu  and  Waters ,  1997]).  We  associate  the  atmospher¬ 
ic  component  (&A)  with  radiance  fluctuations  produced  by 
GWs  and  refer  to  it  hereafter  as  the  GW  variance.  For 
UARS  MLS,  the  instrument  noise  cr^  is  frequency  (chan¬ 
nel)  dependent  but  stable  during  the  entire  mission  [Lau  et 
al ,  1996].  Thus  one  way  to  estimate  the  instrument  noise 
is  to  average  radiance  variances  at  latitudes  (usually  near 
the  tropics)  where  there  is  little  resolved  GW  activity 
(ga  ~  0),  so  that  a2  «  In  practice,  we  estimate  the 
instrument  noise  from  the  minimum  variance  of  monthly 
zonal-mean  averages  (evaluated  globally  within  equispaced 
5°  latitude  bins).  The  minimum  variances  from  about  36 
months  of  such  averages  are  themselves  then  averaged 
together  to  yield  a  mean  estimate  of  the  noise  variance  g2n 
for  each  MLS  channel.  These  noise  variances  (see  Table  1) 
are  within  20%  of  the  calibrated  values  [Wu  and  Waters , 
1997].  GW  variances  then  follow  by  subtracting  the  noise 
variance:  ga  =  b2  —  g2n. 

[n]  The  ga  estimates  are  available  at  the  8  MLS 
channel  saturation  altitudes  (~28,  33,  38,  43,  48,  53,  61, 
and  80  km).  Only  a  small  portion  of  the  atmospheric  GW 
field  is  resolved  by  these  MLS  limb  measurements,  an 


effect  that  must  be  carefully  considered  when  using  ga 
data  to  infer  information  on  fundamental  gravity  wave 
properties  [. Alexander ,  1998;  McLandress  et  al ,  2000]. 
This  “observational  filtering”  is  particularly  acute  for  the 
six-point  variances  we  consider  here  (see,  e.g.,  Figure  1  of 
McLandress  et  al  [2000]).  We  provide  a  thorough  dis¬ 
cussion  of  these  issues  in  section  4. 

2.4.  Further  Improvements 

[is]  In  this  study  we  improved  the  ga  calculations  in 
several  ways.  Specifically,  we  improved  the  quality  of  the 
channel  1  radiance  variances  (which  saturate  at  the  lowest 
altitude),  improved  the  signal-to-noise  ratios  (SNR)  of  all 
the  channel  variances,  and  increased  the  data  volume  that 
went  into  the  mean  variance  statistics. 

[19]  We  found  that  the  tangent  height  zt  must  be  below 
~14  km  in  order  for  channel  1/15  radiances  to  satisfy  the 
saturation  condition.  This  cutoff  limits  the  number  of 
individual  radiance  values  for  the  variance  calculations  in 
equation  (la)  during  normal  limb-scan  operation  to  n  <  4. 
Thus  four-point  limb-scan  variances  replace  the  previous 
six-point  limb-scan  variances  used  by  Wu  and  Waters 
[1996b]  and  give  more  reliable  results  at  lower  altitudes. 
For  limb-track  data,  we  still  compute  six-point  variances 
(zt  ~  18  km),  but  with  the  following  improvements. 

[20]  To  improve  the  SNR  of  ga,  we  combined  radian¬ 
ces  from  a  pair  of  channels  that  are  symmetric  about  the 
line  center  before  performing  the  detrending  operation  in 
equation  (la).  Since  the  symmetric  channels  have  similar 
temperature  weighting  functions  and  mean  noise  levels 
[Wu  and  Waters ,  1997],  the  variance  of  the  combined 
channel  radiances  yields  a  noise  variance  about  half  that 
from  the  single,  channel  radiances  alone,  namely, 
^combined  =  +  “f  •  This  represents  a  factor  of  2  improve¬ 
ment  in  the  variance  SNR.  It  is  important  to  combine  the 
radiances  before  averaging,  otherwise  the  averaged  var¬ 
iances  from  the  two  paired  channels  would  only  improve 
the  SNR  by  a  factor  ofs/l  instead  of  2. 

[21]  Another  improvement  in  the  GW  variance  analysis  is 
to  make  a  full  use  of  all  saturated  radiances.  Instead  of 
cutting  off  the  limb  radiances  at  the  same  tangent  height 
for  all  the  channels,  we  make  the  cutoff  height  channel 
dependent.  Because  radiances  from  different  channels 
saturate  at  different  heights,  higher  cutoff  heights  increase 
the  number  of  samples  p  in  equation  (lb),  reducing  statis¬ 
tical  uncertainties. 

[22]  Finally,  the  new  GW  variance  data  have  more 
accurate  geographical  registration  of  the  atmospheric 
volume  measured  by  each  channel.  Because  saturated 
MLS  radiances  measure  the  atmospheric  volume  on  the 
nearside  of  the  tangent  point  (see,  e.g.,  Figure  2  of  Jiang  et 


Figure  1.  (a)  Variances  from  MLS  channels  3/13  limb-track  measurements  (z  ~  38  km)  during  north-looking  descending 
(ND)  orbit  segments  in  the  Northern  Hemisphere  winter  stratosphere.  These  data  are  averaged  on  a  10°  x  4°  longitude- 
latitude  grid  during  winter  months  (December  through  February)  of  1994-1997.  The  dotted  contour  lines  are  UKMO 
assimilated  mean  stratospheric  winds  at  pressure  altitudes  between  ^6.8  and  4.6  hPa,  averaged  on  the  same  10°,  4° 
longitude-latitude  grid  during  the  same  MLS  limb-track  days  and  interpolated  onto  the  same  limb-track  orbits.  The  color 
bar  scale  for  radiance  variance  is  linear,  (b)  Digital  elevations  of  northern  hemispheric  landmasses,  plotted  using  the 
“ETOP05”  5-min  gridded  digital  elevation  model  issued  by  the  National  Geophysical  Data  Center.  Locations  and  names 
of  mountain  ranges  were  derived  from  various  sources,  the  main  being  National  Geographic  Society  [1995].  The  elevation 
color  bar  scale  is  linear,  (c)  As  for  Figure  la,  but  plotting  variances  on  north-looking  ascending  orbits  (NA). 
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al  [2003]),  the  geographical  location  of  this  atmospheric 
volume  is  displaced  from  the  tangent  point,  in  contrast 
to  the  atmospheric  volume  of  unsaturated  limb  radiances. 
This  correction  was  not  considered  in  previous  MLS 
GW  variance  calculations.  We  found  by  studying  moun¬ 
tain  waves  and  comparing  them  to  model  predictions  in 
this  work  that  inferred  distributions  of  wave  activity 
prove  particularly  sensitive  to  any  such  geolocation 
errors  in  &A. 

3.  Enhanced  MLS  Variances  Over  Mountains 

[23]  To  investigate  possible  variance  enhancements  over 
Northern  Hemisphere  topography,  we  focus  initially  on  the 
limb-track  data,  which  have  higher  along-track  resolution 
and  are  more  easily  analyzed  for  their  GW  content  than 
limb-scan  data  [McLandress  et  al ,  2000].  Figure  la  plots  a 
hemispheric  map  of  mean  six-point  limb  track  radiance 
variances  from  channels  3  and  13  (altitude  ^38  km)  during 
northern  winter  (December-February)  for  the  years  1994- 
1997,  as  acquired  during  north  looking  (N)  yaw  cycles  and 
descending  (D)  orbit  segments.  We  refer  to  these  hereafter 
as  “ND  variances,”  ( ®a)nd •  NA  variances  (a a)na  derived 
from  north  looking  (N)  ascending  (A)  orbits  are  plotted  in 
Figure  lc.  Figure  lb  plots  a  high-resolution  map  of  the 
topographic  elevations  of  landmasses  in  the  Northern  Hemi¬ 
sphere,  with  major  features  relevant  to  potential  mountain 
wave  forcing  highlighted.  This  includes  major  continental 
mountain  ranges,  mountainous  islands,  and  the  Greenland 
mountain-ice  complex.  The  major  focus  is  topography 
poleward  of  30°N,  since  MLS  variances  equatorward  of 
30°  in  Figures  la  and  lc  show  no  significant  variance 
enhancements. 

3.1.  Sensitivity  of  Variances  to  MLS  View  Direction 

[24]  MLS  can  view  low-latitude  regions  of  the  atmo¬ 

sphere  in  any  one  of  four  different  viewing  geometries, 
which  vary  with  the  north-south  yaw  cycle  of  UARS  and 
whether  the  spacecraft  is  on  the  ascending  (south  to  north) 
or  descending  (north  to  south)  portion  of  its  orbital  track. 
Only  north-viewing  yaw  cycles  yield  near-hemispheric 
coverage  of  the  Northern  Hemisphere,  which  limits  us  to 
the  two  ND  and  NA  views  introduced  in  Figures  la  and  lc, 
respectively.  ND  and  NA  variances  differ  significantly. 
Most  obvious  is  an  overall  decrease  in  intensity  and  a  much 
“spottier”  global  distribution  of  (ga)na  compared  to 
(a a) nd •  et  al-  [2003]  found  a  similar  result  on  plotting 

MLS  variance  maps  over  the  southern  Andes  and  western 
USA.  The  visual  correlations  with  topography  (Figure  lb) 
are  much  more  striking  for  (g2a ) ND  in  Figure  la  than  for 

(°a)na  in  Fi§ure  lc- 

[25]  This  sensitivity  to  view  direction  was  studied  and 
quantified  somewhat  by  McLandress  et  al  [2000].  Why  it 
arises  physically  is  depicted  schematically  in  Figure  2,  for 
the  particular  case  of  a  stationary  mountain  wave  in  the 
northern  winter  midlatitude  stratosphere.  Predominantly 
eastward  winds  generate  mountain  waves  propagating  west¬ 
ward  with  respect  to  this  flow,  yielding  waves  with  three- 
dimensional  (3-D)  phase  fronts  aligned  broadly  as  depicted 
in  Figure  2.  Various  MLS  limb  view  directions,  labeled  1 
though  5,  are  depicted  above  this  wave.  Each  limb  view 
intercepts  the  3-D  wave  structure  at  the  indicated  slant  angle 


Figure  2.  Schematic  3-D  depiction  of  mountain  waves 
forced  by  horizontal  flow  (U,  V)  over  an  idealized  mountain 
ridge.  The  wind  profile  (depicted  with  hatched  vectors) 
remains  largely  eastward  through  the  stratosphere  (U  >  0), 
supporting  stratospheric  mountain  waves  with  horizontal 
wave  numbers  directed  largely  westward  (k  <  0,  cp  ~  180°) 
with  three-dimensional  phase  fronts  oriented  broadly  as 
depicted.  Five  potential  MLS  slant  viewing  geometries 
(slant  angle  (3),  spanning  180°  in  viewing  azimuth  aT  and 
labeled  1  through  5,  are  shown  above  this  mountain  wave 
and  are  envisaged  as  passing  through  its  three-dimensional 
structure  while  acquiring  saturated  radiances  from  this 
volume  of  atmosphere  directly  above  the  ridge.  Viewing 
geometry  2  is  most  favorable  for  resolving  this  wave 
perturbation  in  saturated  along-track  radiances,  since  it 
views  quasi-parallel  to  the  wave’s  phase  fronts  and  projects 
a  long  horizontal  wavelength  along  and  toward  the  LOS, 
yielding  less  along-path  smearing  of  the  wave  oscillation. 
Conversely,  view  5  views  the  limb  across  the  wave’s  phase 
fronts  and  so  will  tend  to  smear  out  the  wave  in  each 
instantaneous  limb  radiance  measurement. 


(3  and  varying  azimuth  angles  aT  in  acquiring  saturated 
radiances  from  this  specific  volume  of  the  atmosphere.  MLS 
view  directions  1-3  in  Figure  2  are  the  most  favorable  for 
resolving  this  wave,  since  they  intercept  the  wave  roughly 
parallel  to  its  sloping  3-D  wave  fronts,  whereas  views  4-5 
sample  it  at  progressively  more  anti-parallel  alignments 
with  respect  to  wave  phase,  which  smears  the  wave’s 
signature  out  along  the  limb.  For  north  yaw  cycles,  MLS 
views  toward  the  east  during  descending  orbits,  yielding 
viewing  geometries  somewhat  like  those  labeled  1  and  2  in 
Figure  2,  whereas  ascending  orbits  produce  westward  MLS 
view  directions  and  unfavorable  viewing  geometries  more 
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like  those  labeled  4-5  in  Figure  2.  For  further  illustrations 
of  these  effects,  see  Jiang  et  al.  [2003]. 

[26]  For  these  reasons,  we  focus  primarily  on  (v2a)nd  in 
our  initial  investigation  of  possible  mountain  wave  signals 
in  MLS  northern  hemispheric  radiances.  However,  this 
simple  ND/NA  discrimination  is  just  the  first  step  in 
unraveling  the  complex  ways  in  which  MLS  viewing  angles 
influence  the  mountain  wave  content  of  the  dA  data.  We 
defer  a  thorough  treatment  of  these  issues  to  section  4.2  as 
part  of  our  attempt  to  model  the  mountain  wave  content  of 
both  the  {v2a)nd  and  (<j2a ) NA  data  globally. 

3.2.  Sensitivity  to  Background  Wind  Speed 

[27]  Dotted  contours  in  Figures  la  and  lc  show  mean 
horizontal  wind  speeds  from  the  United  Kingdom  Meteo¬ 
rological  Office  (UKMO)  assimilations  [Swinbank  and 
O’Neill ,  1994],  averaged  for  the  specific  MLS  observation 
days  that  contributed  data  to  these  variances.  Figure  la 
shows  a  strong  correlation  between  enhanced  (Pa)nd  and 
strong  stratospheric  horizontal  wind  speeds  U  that  has  been 
noted  in  previous  studies  [e.g.,  Wu  and  Waters ,  1996a, 
1996b,  1997;  Jiang  and  Wu ,  2001].  Regions  in  Figure  la 
where  U  is  less  than  ^10  m  s-1  correspond  to  very  small 
(d2A)ND  (mainly  blue  contours  where  radiance  variances  are 
<0.01  K2).  Conversely,  regions  of  largest  (cr2 )  ND  all  occur 
in  midlatitude  to  high-latitude  regions  where  mean  strato¬ 
spheric  winds  U  are  ~25  ms-1  or  greater.  These  {®2a)nd 
enhancements  follow  to  some  extent  the  zonal  modulation 
of  the  flow  by  the  climatological  quasi- stationary  wave 
number- 1  pattern  that  emerges  in  these  UKMO  wind 
averages,  such  that  strong  (weak)  midlatitude  mean  winds 
over  Eurasia  (North  America)  correspond  to  similarly  strong 
(weak)  (a2A )  ND  over  each  continent. 

[28]  Alexander  [1998]  showed  that  the  primary  source  of 
this  U-  (cf 2 )  VD  correlation  is  the  refraction  of  GW  vertical 
wavelengths  by  background  winds,  coupled  with  the  deep 
vertical  width  of  MLS  weighting  functions  [Wu  and  Waters , 
1997]  that  do  not  allow  MLS  to  resolve  GWs  with  vertical 
wavelengths  \z  <  10  km  or  so.  For  the  specific  case  of 
stationary  hydrostatic  mountain  waves,  the  vertical  wave¬ 
length  \z  is  given  by 


where  N  is  the  Brunt- Vaisala  frequency  and  0  is  the  angle 
between  the  horizontal  wind  U  and  the  wave’s  horizontal 
wave  number  vector  Kh  (see,  e.g.,  Figure  2).  At  latitudes 
~40°-60°N  over  central  Eurasia,  U  ~  40  to  50  ms-1  from 
Figure  la.  If  we  assume  for  now  that  some  of  these 
mountain  waves  are  aligned  parallel  to  the  mean  flow  (0  ~ 
0°  or  180°),  then  N  ~  0.02  rad  s-1  so  that  (3)  yields  Xz  ~ 
13-16  km.  Thus  MLS  should  resolve  these  waves, 
consistent  with  the  variance  enhancements  observed  here 
in  Figure  la.  Conversely,  a  similar  calculation  for  mountain 
waves  over  the  Himalayas  at  ^30°N,  where  U  <  10  ms-1, 
yields  \z  <  3  km,  a  vertical  wavelength  too  short  for  any 
MLS  channel  to  resolve. 

3.3.  Enhanced  (oa)nd  Over  Mountainous  Regions 

[29]  The  analysis  in  sections  3.1  and  3.2  has  introduced 
some  important  ways  in  which  mountain  waves  may  and 


may  not  manifest  in  MLS  limb  track  radiances.  With  these 
issues  in  mind,  we  now  compare  Figures  la  and  lb  and 
identify  four  geographical  “focus  regions”  with  significant 
underlying  topography  where  enhanced  stratospheric 
(Pa)nd  ^eve^s  are  evident.  We  choose  these  particular 
regions  on  the  basis  of  previous  observational  precedents 
for  supposing  that  these  variance  enhancements  are  good 
potential  candidates  for  a  mountain  wave  explanation. 

[30]  1.  The  first  focus  region  is  the  Kjonas  Mountains  of 
Scandinavia  (~20°E,  65 °N):  Clear  {g2a)nd  bursts  are  evi¬ 
dent  in  Figure  la  along  the  northern  and  southern  ends  of 
these  mountains.  Regional  observations  and  modeling  stud¬ 
ies  have  shown  that  eastward  flow  from  the  Arctic  Ocean 
across  these  mountains  during  winter  can  generate  strato¬ 
spheric  mountain  waves  with  long  horizontal  and  vertical 
wavelengths  [e.g.,  Dornbrack  et  al ,  2002]  that  should 
theoretically  be  detectable  in  satellite  limb  radiances.  To 
study  MLS  variances  here  further,  we  define  a  latitude- 
longitude  box  of  55°-70°N,  5°-40°E. 

[31]  2.  The  second  focus  region  is  Central  and  Western 
Eurasia  (~100°E,  50°N):  Large  geographically  extended 
(52)Nd  values  occur  in  this  region  where  a  broad  complex 
of  interconnected  mountain  ranges  exist:  the  Altai,  Hangay 
and  Sayan  Mountains  to  the  west,  the  Yablonovyy,  Khin- 
gan,  and  Stanavoy  Ranges  to  the  east.  Eckermann  and 
Preusse  [1999]  noted  stratospheric  variance  enhancements 
in  a  similar  location  during  CRISTA  infrared  stratospheric 
limb  observations  in  November  1994  and  used  the  MWFM 
(see  section  4.1)  to  identify  them  as  stratospheric  mountain 
waves  forced  by  flow  across  these  ranges.  We  define  a 
latitude-longitude  box  in  the  range  45°-65°N,  80°-130°E 
to  study  variances  in  this  region  further. 

[32]  3.  The  third  focus  region  is  the  Queen  Elizabeth 
Islands  (~270°E,  70°N):  These  Canadian  Arctic  islands 
include  regions  of  significant  topography,  such  as  Baffin 
and  Ellesmere  Islands,  the  latter  slightly  poleward  of  the 
most  northern  observing  latitudes  for  saturated  MLS  radi¬ 
ances.  Accumulated  high-resolution  lidar  and  balloon  pro¬ 
filing  of  the  winter  stratosphere  over  stations  in  this  region 
has  revealed  gravity  waves  whose  temperature  variances 
vary  in  intensity  with  background  winds  in  ways  consistent 
with  quasi-stationary  phase  speeds,  and  thus  wave  energy 
here  has  been  attributed  mostly  to  locally  generated  moun¬ 
tain  waves  [e.g.,  Whiteway  and  Duck ,  1996,  1999;  Duck  et 
al ,  1998;  Duck  and  Whiteway ,  2000].  For  further  study  of 
(d2)ND  °ver  these  islands,  we  define  a  latitude-longitude 
box  in  the  range  65°-75°N,  80°-120°W. 

[33]  4.  The  fourth  focus  region  is  Southern  Greenland 
(^315°E,  65°N):  In  addition  to  the  large  ice  cap  that 
dominates  Greenland’s  topographic  relief,  significant  nuna- 
taks  (mountains  that  protrude  through  glacial  ice)  occur  all 
along  the  east  coast  to  its  southernmost  tip.  An  ER-2  flight 
over  southern  Greenland  on  6  January  1992  recorded 
enhanced  velocity  fluctuations  in  the  stratosphere,  which 
Leutbecher  and  Volkert  [2000]  simulated  with  a  mesoscale 
model  and  attributed  to  mountain  waves  forced  by  flow 
over  southern  Greenland’s  topography.  To  study  (a in 
this  region  further,  we  define  we  define  a  latitude-longitude 
box  in  the  range  55°-70°N,  20°-60°W. 

[34]  In  studying  (cr a)nd  variability  further  over  our  four 
mountainous  focus  regions  defined  above,  we  fold  in  limb- 
scan  data  in  addition  to  limb-track  data.  UARS  MLS 
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acquired  ^6  years  of  limb-scan  data  (1991-1997)  in 
addition  to  the  ^3  years  of  limb-track  data  (1994- 
1997).  The  differences  between  limb-scan  and  limb-track 
data  are  mainly  in  viewing  geometry  and  sampling 
periods  [McLandress  et  al .,  2000;  Jiang  et  al .,  2003]. 
The  different  viewing  geometries  give  rise  to  systematic 
differences  in  the  radiance  variances  computed  from  each 
data  type.  For  example,  four-point  and  six-point  limb-scan 
variances  are  generally  much  larger  than  corresponding 
six-point  limb-track  variances,  which  McLandress  et  al 
[2000]  attributed  to  aliasing  of  longer  wavelength  GWs 
into  limb-scan  data  due  to  a  varying  LOS  angle,  effects 
we  will  investigate  further  in  subsequent  studies. 

3.3.1.  Seasonal  and  Interannual  Variations 

[35]  Figure  3  plots  time  series  of  (®2a)nd  for  channels 
1/15,  3/13  and  5  (pressure  altitudes  z  ^  28  km,  ~38  km 
and  ^48  km,  respectively)  over  our  four  mountainous 
areas  of  interest  just  defined  (plotted  on  the  inset  map  in 
Figure  3f):  Limb-scan  variances  are  shown  in  Figures  3 a- 3c, 
limb-track  variances  in  Figures  3d-3f.  At  all  four  focus 
regions  we  find  a  strong  reproducible  annual  cycle,  peaking 
in  winter. 

[36]  Rocket  measurements  of  stratospheric  GW  tempera¬ 
ture  variances  at  middle  and  high  northern  latitudes  also 
show  annual  variations  peaking  in  winter  [e.g.,  Hirota  and 
Niki ,  1985;  Eckermann  et  al ,  1995],  as  do  GPS/MET 
satellite  measurements  [Tsuda  et  al ,  2000].  At  high  lati¬ 
tudes,  Eckermann  [1995]  showed  that  a  seasonal  variation 
in  atmospheric  densities  at  a  given  geometric  altitude  could 
reproduce  most  of  the  observed  annual  cycle  in  the  rocket 
data.  As  noted  by  Alexander  [1998],  a  similar  model  cannot 
explain  the  annual  variations  of  (cbj )  ND  in  Figure  3,  since 
MLS  channel  altitudes  are  registered  at  pressure  heights 
rather  than  geometric  heights,  which  largely  factors  out  the 
seasonal  density  effect.  Using  a  ray  model,  Alexander 
[1998]  attributed  much  of  the  annual  variability  in  MLS 
variances  to  seasonal  variations  in  stratospheric  winds  that 
allowed  more  waves  in  winter  to  attain  long  vertical 
wavelength  and  become  visible  to  MLS,  in  the  manner 
outlined  in  section  3.2. 

[37]  In  their  study  of  MLS  variances  over  the  southern 
Andes,  Jiang  et  al  [2002]  also  found  an  annual  variation, 
peaking  in  winter,  which  they  attributed  via  modeling 
studies  to  increased  mountain  wave  activity  in  winter.  The 
variance  enhancements  observed  at  our  four  northern  hemi¬ 
sphere  focus  regions  also  occur  during  winter  and  persist 
into  early  spring,  abating  rapidly  in  April.  Very  similar 
winter-to-spring  variations  have  been  observed  in  strato¬ 
spheric  GWs  profiled  at  stations  in  and  around  the  Queen 
Elizabeth  Islands,  where  enhanced  winter  variances  corre¬ 
late  with  strong  eastward  vortex  winds  [e.g.,  Whiteway , 
1999;  Whiteway  and  Duck ,  1999].  Sudden  reductions  in 
early  spring  occur  with  final  breakdown  of  the  winter 


vortex,  which  yields  weaker  winds  and  reduced  transmis¬ 
sion  of  locally  generated  mountain  waves  into  the  strato¬ 
sphere  [ Duck  and  Whiteway ,  2000].  Much  reduced 
variances  in  summer  are  consistent  with  westward  flow 
which  prevents  mountain  wave  penetration  into  the  mid¬ 
dle  stratosphere  via  critical-level  filtering.  We  also  note 
considerable  site-to-site  and  interannual  variability  of  the 
wintertime  peaks  in  Figure  3,  which  is  consistent  with  the 
interannual  variability  of  the  winter  vortex  which  gives 
rise  to  considerable  interannual  variations  in  the  intensity 
of  Arctic  mountain  waves  entering  the  stratosphere 
[Dornbrack  and  Leutbecher ,  2001].  However,  larger 
limb-scan  variances  after  1995  are  an  observational  arti¬ 
fact  corresponding  to  changes  from  “normal  scan”  mode 
(top  to  bottom)  to  “reverse  scan”  mode  (bottom  to  top), 
the  latter  yielding  generally  larger  (^a)nd  variances- 

[38]  In  summary,  time  series  of  ( ®a)nd  Figure  3  are 
consistent  with  increases  in  resolved  stratospheric  mountain 
waves  over  the  four  focus  regions  during  winter.  This 
interpretation  is  not  definitive,  however,  since  complex 
wind-modulated  visibility  effects  outlined  in  section  3.2 
yield  a  similarly  phased  annual  variation  [„ Alexander ,  1998]. 

3.3.2.  ND  Versus  NA  Variances 

[39]  As  noted  in  section  3.1,  mountain  waves  in  winter 
should  manifest  more  strongly  in  {ga)nd  relative  to  {d2A)NA. 
Figures  la  and  lc  showed  that  climatologically,  mean  ND 
limb-track  variances  in  northern  winter  were  indeed  signif¬ 
icantly  larger  than  the  NA  variances  over  most  mountainous 
regions.  To  study  this  in  more  depth,  Figure  4  shows 
scatterplots  of  daily  mean  ND  versus  NA  limb-scan  var¬ 
iances  over  our  four  mountainous  regions  of  interest  during 
winter  months  from  1991  to  1997.  When  variances  are 
strongly  enhanced  (i.e.,  M).l  K2),  we  note  a  much  greater 
likelihood  that  this  enhanced  variance  is  observed  during 
ND  views,  whereas  NA  views  of  the  same  region  yield 
moderate  or  small  variances.  This  bias  is  consistent  with 
what  we  would  expect  if  MLS  fluctuations  over  these 
regions  were  enhanced  by  viewing  large-amplitude  moun¬ 
tain  waves  propagating  southwestward  in  a  mainly  eastward 
winter  stratospheric  flow  (as  depicted  in  Figure  2). 

3.3.3.  Height  Variations 

[40]  Figure  5  plots  mean  altitude  profiles  of  limb-scan 
(d2A)ND  normalized  by  mean  radiance  brightness  temper¬ 
atures,  derived  from  the  6  years  of  measurements  in  all 
channels  above  the  four  mountainous  focus  regions  during 
winter.  These  regional  profiles  resemble  zonal-mean  pro¬ 
files  at  high  latitudes  reported  in  previous  studies  [Wu  and 
Waters ,  1996b;  Wu  and  Jiang ,  2002]:  specifically,  exponen¬ 
tial  growth  up  to  an  altitude  ^50  km,  then  an  apparent 
“saturation”  of  the  variance  profile  at  altitudes  above  50  km. 
Values  here  are  somewhat  smaller  than  those  reported  by 
Wu  and  Waters  [1996b]  because  of  a  four-point  rather  than 
six-point  variance  calculation. 


Figure  3.  Time-series  of  (a-c)  daily  mean  MLS  limb-scan  GW  variances  for  normal-scan  period  (October  1991  to 
December  1994),  reverse-scan  period  (January  1995  to  May  1997),  and  (d-f)  daily  mean  MLS  limb-track  GW  variances 
from  December  1994  to  May  1997.  The  variances  are  from  UARS  MLS  north-looking  descending  orbits  at  three  different 
altitudes  of  28  km  (channels  1/15),  38  km  (channels  3/13),  and  48  km  (channels  5)  over  four  selected  focus  regions  shown 
on  the  inset  map  to  the  lower  left  of  Figure  3f:  southern  Greenland  (55°N-70°N;  20°-60°W;  blue);  central  Eurasia  (45°N- 
65°N;  80°-130°E;  red),  Scandinavia  (55°N-70°N;  5°-40°E;  green);  Queen  Elizabeth  Islands  (65°N-75°N;  80°-120°W; 
black).  A  3 -day  running  mean  was  applied  to  the  daily  averages  to  smooth  out  the  large  fluctuations  in  the  daily  variances. 
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Greenland  Scandinavia 


Figure  4.  (a-d)  Scatterplot  of  daily  mean  variances  at  the 
four  focus  regions,  derived  from  limb-scan  radiance 
measurements  (1991-1997)  from  channels  3  and  13  (z  ~ 
38  km)  on  north-looking  descending  (ND)  and  ascending 
(NA)  orbits  for  winter  months  (December- February).  The 
variances  are  computed  in  the  same  way  and  for  the  same 
four  regions  as  in  Figure  3. 


[41]  Wu  and  Waters  [1996b,  1997]  interpreted  variations 
below  50  km  in  terms  of  the  ~exp(J  dz/HE)  growth  in 
wave-induced  temperature  variances  with  height  [Fritts  and 
VanZandt ,  1993;  Eckermann,  1995].  When  He  ~  Hp  «  7  km 
(where  Hp  is  density  scale  height),  all  the  waves  are  non¬ 
dissipating:  Growth  at  this  rate  is  shown  with  a  grey  curve 
in  Figure  5.  “Energy  scale  heights”  HE  for  the  MLS 
variances  below  ^50  km  in  Figure  5  are  ~12  km,  sugges¬ 
tive  of  some  limited  dissipation  of  wave  energies  with 
height. 

[42]  Above  ~50  km,  variances  in  Figure  5  become  almost 
constant  with  height  (HE  — ►  00).  Two  interpretations  of  this 
observation  have  been  proposed.  Wu  and  Waters  [1996b, 
1997]  interpreted  this  as  a  transition  to  greater  GW  dissi¬ 
pation,  leading  to  saturation  limits  at  high  altitudes  indica¬ 
tive  of  strong  wave  breaking.  Wu  and  Jiang  [2002] 
reiterated  this  interpretation  when  studying  variances  over 
Antarctica,  while  Wu  [2001]  proposed  a  similar  explanation 
for  vertical  growth  tendencies  observed  in  along-track  wave 
number  spectra  of  MLS  radiances.  Alexander  [1998],  how¬ 
ever,  used  a  global  GW  model  to  argue  that  wind-modulated 
MLS  visibility  effects  in  section  3.2  produced  most  of  the 
apparent  “saturation”  in  variances  above  50  km,  and  that 
profiles  like  those  in  Figure  5  gave  little  indication  of  where 
the  strongest  wave  breaking  was  occurring. 

[43]  We  leave  these  interpretation  issues  for  future  dedi¬ 
cated  modeling  studies.  Here,  we  point  out  that  the  vertical 


variations  of  MLS  variances  over  these  four  mountainous 
locations  are  generally  consistent  with  earlier  zonal  mean 
results  and  indicate  that  the  GW  activity  here  propagates 
through  the  stratosphere  and  into  the  mesosphere.  This  is 
consistent  with  a  mountain  wave  interpretation,  since  east¬ 
ward  winter  stratospheric  flow  persists  climatologically  into 
the  lower  mesosphere  and  so  should  often  allow  mountain 
waves  to  propagate  through  the  full  depth  of  the  winter 
middle  atmosphere  from  0  to  50  km.  The  strong  altitude 
correlations  between  winter  variance  enhancements  in 
Figure  3  further  support  such  conclusions. 

4.  MWFM  Modeling  of  Mountain  Waves 
in  MLS  Variances 

[44]  Data  analysis  in  section  3  has  provided  evidence  that 
mountain  waves  may  be  enhancing  (cf^ )  ND  variances  over 
topography.  To  investigate  these  issues  in  depth,  we  develop 
here  a  global  modeling  strategy  for  estimating  the  mountain 
wave  content  of  MLS  radiances  using  the  Mountain  Wave 
Forecast  Model  (MWFM).  We  begin  by  briefly  introducing 
the  MWFM,  develop  in  detail  a  critical  “MLS  visibility 
function,”  which  specifies  the  waves  that  MLS  can  and 
cannot  see  in  orbit,  and  then  apply  it  to  some  simple 
examples.  Finally  we  outline  details  of  our  MWFM  simu¬ 
lations  with  this  visibility  term  included  that  are  targeted  to 
the  MLS  observations  presented  in  section  3. 

4.1.  Mountain  Wave  Forecast  Model  (MWFM) 

[45]  Bacmeister  et  al.  [1994]  documents  the  first  version 
of  this  model  (MWFM-1).  Here,  as  in  the  study  by  Jiang  et 
al  [2002],  we  use  version  2  (MWFM-2).  The  MWFM  uses  a 
set  of  diagnosed  quasi-two  dimensional  (2-D)  ridge  func¬ 
tions  that  fit  the  dominant  features  in  the  Earth’s  topography 
relevant  for  mountain  wave  generation  [Bacmeister  et  al , 


Figure  5.  Vertical  growth  of  variances  over  southern 
Greenland,  central  Eurasia,  Scandinavia,  and  Queen  Eliza¬ 
beth  Islands  for  Arctic  winters  (1991-1997)  as  observed  in 
limb  scans  by  MLS  on  north-looking  descending  (ND) 
orbits.  The  variances  are  normalized  by  the  mean-square 
radiance  brightness  temperature.  Instrument  noise  was 
removed  from  these  variances.  The  grey  curve  shows 
sample  exp  (J  dz/H p)  growth  with  H=  7  km. 
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1994].  Each  ridge  function  has  a  range  of  properties:  The 
most  important  are  its  geographical  location,  cross-ridge 
width  L ,  peak  height  h ,  base  altitude  zb  and  horizontal 
orientation  cp LONG  of  its  long  (along-ridge)  axis.  Also  im¬ 
portant  for  MWFM-2  calculations  is  a  normalized  “ridge 
quality”  parameter  g  (0  <  q  <  1),  which  defines  how  closely 
the  original  topography  was  approximated  by  fitting  a  2-D 
ridge  function:  q  ~  1  indicates  highly  2-D  ridge-like 
topography,  whereas  g  C  1  usually  indicates  a  more 
symmetric  three-dimensional  (3-D)  obstacle. 

[46]  The  model  ingests  forecast  or  analyzed  wind  fields 

and  “blows”  the  local  winds  at  the  base  altitude  zb  of  the 
ridge  across  it  to  generate  a  source-level  mountain  wave 
field.  MWFM-2  models  forced  mountain  wave  fields  using 
a  nonhydrostatic  three-dimensional  spatial  ray-tracing  for¬ 
mulation  with  Coriolis  and  stratification  terms  included 
[Marks  and  Eckermann ,  1995],  constrained  by  a  uniformly 
stationary  ground-based  phase  speed  for  each  ray.  A  series 
of  rays  are  launched  from  the  parent  ridge  with  a  given 
horizontal  wave  number  vector  (Kb ) .  directed  at  fy  different 
equispaced  azimuth  angles  cpz  .  This  set  of  azimuth  angles 
span  a  180°  range  g>OR  —  90°  <  <  g)OR  +  90°,  where  g)OR 

is  the  azimuth  orthogonal  to  the  long  axis  of  the  ridge 
(i p0R  =  long  =*=  90°)  and  i  is  an  integer  that  labels  each  ray 

O'  —  1  . . .  I). 

[47]  Rays  aligned  parallel  to  g)OR  are  forced  with  the 

largest  peak  vertical-displacement  amplitude  based  on  the 
peak  ridge  height  h.  The  ridge  quality  q  is  used  to  scale  wave 
amplitudes  down  from  this  value  as  a  function  of  displace¬ 
ment  of  the  ray  azimuth  from  \g>t  —  For  q  ~  0 
(quasi-3-D  obstacle),  the  amplitude  reduction  as  a  function 
of  |cpz-  —  is  small,  whereupon  the  multiray  method  often 
yields  three-dimensional  “ship  wave”  patterns  radiating 
from  the  obstacle  as  rays  at  various  azimuths  radiate 
obliquely  away  from  the  source  with  significant  amplitudes 
[see,  e.g.,  Gjevik  and  Marthinsen,  1978;  Broutman  et  al, 
2001].  For  q  ~  1  (2-D  ridge-like  topography),  initial  wave 
amplitudes  decay  rapidly  with  increasing  |cp-  —  | ,  and 

thus  the  ray  method  yields  significant  energy  only  for  those 
few  rays  aligned  close  to  g)OR,  yielding  a  plane  wave  pattern 
with  phase  lines  quasi-parallel  to  the  long  ridge  axis,  much 
as  is  observed  for  flows  over  two-dimensional  ridges  (see 
Figure  2).  Once  assigned,  these  wave  amplitudes  vary  along 
each  ray’s  group  propagation  path  according  to  conservation 
of  the  vertical  flux  of  wave  action  density,  subject  to 
dissipation  by  dynamical  and  convective  wave  breaking 
thresholds  [Marks  and  Eckermann ,  1995].  Wave  breaking 
is  parameterized  using  a  linear  saturation  hypothesis,  with 
wave  action  densities  scaled  back  to  saturation  thresholds 
where  breaking  occurs.  At  any  altitude  the  wave  action 
density  can  be  converted  to  a  more  usual  wave  amplitude 
measure  [Marks  and  Eckermann ,  1995].  In  these  experi¬ 
ments  we  convert  to  peak  temperature  amplitude  7),  which 
relates  most  closely  to  the  radiance  fluctuations  measured  by 
MLS. 

[48]  The  total  number  of  rays  launched  from  each  ridge 
/  =  fyfy,  where  IK  is  the  total  number  of  horizontal 
wave  numbers  assigned  to  each  ridge  and  is  the 
total  number  of  propagation  azimuths  per  horizontal  wave 
number.  In  the  experiments  reported  here,  |(^).|  = 
(Kh)i~  1  -5/, //-(./,  —  1  . . .  /A-.  L  is  cross  ridge  width),  IK 
=  2  and  fy  =  18,  yielding  a  total  of  36  rays  per  ridge.  Since 


we  ignore  horizontal  gradient  terms  in  our  ray  formulation, 
these  source-level  horizontal  wave  numbers  and  azimuths 
remain  constant  along  each  ray  trajectory. 

[49]  In  the  MWFM-2  we  use  the  sign  convention  of 
Eckermann  [1992].  We  define  m  <  0  for  upward  propa¬ 
gating  mountain  waves.  For  the  horizontal  wave  number 
vector  (Eh)j=  (ktj  /*■)  =  (Kb)i(cos  <pz-,  sincpz-),  we  define  the 
absolute  wave  number  azimuth  cpz-  such  that  the  compo¬ 
nent  of  (Kh)i  in  the  direction  of  the  background  hori¬ 
zontal  wind  vector  U  =  (£/,  V)  =  C4(cosX;  sin\)  is 
always  negative.  This  always  yields  ( K^j.-U  <  0  and 
thus  positive  intrinsic  frequencies  (ut  =  -(kh).-u  >  0)  ,  so 
that  the  intrinsic  vertical  phase  velocity  //wz  <  0  and 
thus  phase  moves  downward  for  waves  propagating 
energy  upward.  Further,  intrinsic  horizontal  phase  speed 
along  the  horizontal  wind  vector  direction  is  wi/[(Kb)i 
cos(cpz-  —  x)]  =  —Uh,  which  has  the  mountain  wave  prop¬ 
agating  its  phase  in  the  intrinsic  frame  opposite  to,  but  at 
the  same  speed  as  the  wind,  so  that  phase  in  the  ground- 
based  frame  is  always  stationary,  as  required. 

[50]  This  consistent  sign  convention  proves  particularly 
important  for  quantifying  the  precise  sensitivity  of  MLS 
to  MWFM-simulated  mountain  waves,  as  will  now  be 
described. 

4.2.  MLS  Analytical  Filter  Function 
and  Implementation  in  MWFM 

[51]  The  MWFM  generates  mountain  waves  with  a  broad 
range  of  possible  vertical  and  horizontal  wavelengths,  yet, 
as  noted  in  sections  3.1  and  3.2,  MLS  can  only  resolve  a 
small  subset  of  these  waves.  To  enable  meaningful  MWFM- 
MLS  comparisons,  we  must  apply  an  observational  filter  to 
the  MWFM  results  that  accurately  specifies  the  sensitivity 
of  MLS  to  mountain  waves  of  various  horizontal  and 
vertical  wavelengths  aligned  at  various  angles  to  the  MLS 
LOS. 

[52]  In  earlier  MWFM  modeling  of  MLS  variances  over 
the  Andes,  Jiang  et  al  [2002]  adopted  a  simple  filter  that 
retained  waves  with  vertical  wavelengths  \z  >  10  km  and 
horizontal  wavelengths  \b  >  30  km,  similar  to  the  earlier 
modeling  approach  of  Alexander  [1998].  However,  on 
hemispheric  scales,  Wu  and  Waters  [1997]  and  McLandress 
et  al  [2000]  showed  that  there  was  an  important  azimuthal 
sensitivity  to  the  MLS  response  that  changed  as  the  MLS 
viewing  angle  changed  from  the  bottom  (top)  to  the  top 
(bottom)  of  any  given  ascending  (descending)  orbital  seg¬ 
ment.  McLandress  et  al  [2000]  derived  an  analytical 
observational  filter  function  for  MLS  limb  track  observa¬ 
tions,  which  approximated  the  global  three-dimensional 
response  to  GWs.  In  order  to  simulate  how  MWFM- 
simulated  mountain  waves  in  the  NH  winter  stratosphere 
manifest  in  MLS  data,  we  developed  a  fairly  general  MLS 
filter  function  that  builds  upon  this  existing  work,  and 
implemented  it  within  the  MWFM.  Given  its  important  role 
in  model-data  comparisons,  a  frill  description  follows.  We 
do  not  derive  a  similar  filter  function  for  the  limb  scan  data 
here:  As  noted  by  McLandress  et  al  [2000],  additional 
complications  (e.g.,  variable  LOS)  makes  analytical  treat¬ 
ment  of  the  limb  scan  response  considerably  more  difficult. 

[53]  Figure  6  shows  the  horizontal  geometry  of  the  MLS 
observation  for  a  descending  orbital  segment.  Given  our 
focus  on  ND  variances  (cr^)™  we  consider  the  north/ 
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MIS  Viewing  Geometry:  Descending  Orbit 


Figure  6.  Schematic  diagram  of  the  in-orbit  geometry  of 
the  MLS  observation  for  a  north/descending  (ND)  orbit 
segment.  See  text  for  additional  details. 


approach  here  is  to  treat  the  yaw  maneuver  as  a  180° 
rotation  of  the  (X,  Y)  coordinate  axes.  The  advantage  is  that 
this  and  all  subsequent  axis  transformations  are  rotations 
and  so  right-handed  coordinate  axes  persist  throughout.  The 
disadvantage  is  that  the  X  axis  is  rotated  such  that  for 
southern  views,  the  satellite  motion  is  in  the  negative 
direction  along  the  X  axis,  rather  than  in  the  positive 
direction  specified  in  the  original  aT  definition.  The  MLS 
visibility  functions  introduced  below  prove  insensitive  to 
the  absolute  sense  of  quantities  along  the  X  axis,  so  for  our 
purposes  this  side  effect  poses  no  problems. 

[55]  The  MLS  view  along  the  Y  axis  points  downward 
through  the  atmosphere  at  an  angle  (3  with  respect  to  the 
horizontal  plane,  where  by  geometry  (see  Figure  lib  of 
McLandress  et  al  [2000]) 

(3  =*?*-  cos-1  [{a  +  zt)/{a  +  zc)\  (5) 

and  the  radius  of  the  Earth  a  =  6378  km.  For  channels  3  and 
13,  (3  =  —4.53°  and  thus  this  small  angle  allows  us  to 
assume  sin  (3  «  (3  to  a  good  approximation. 

[56]  The  vertical  temperature  weighting  functions  for 
each  MLS  channel  [Wu  and  Waters ,  1996b]  are  approxi¬ 
mated  using  the  Gaussian  fit 


descending  MLS  view,  shown  in  darker  colors  on  the  top- 
right.  MLS  limb-track  radiances  for  a  given  channel  saturate 
at  some  altitude  zc\  For  channels  3  and  13,  zc  =  38  km  and 
the  tangent  height  zt=  18  km  [Wu  and  Waters ,  1997].  The 
portion  of  the  atmosphere  centered  at  zc  contributes  most  to 
the  retrieved  MLS  radiance,  and  so  we  define  this  as  the 
instantaneous  “data  point”  and  the  line  of  data  points 
acquired  following  the  satellite  motion  as  the  “data  track”. 
We  define  (x,  y,  Z)  axes  with  the  origin  at  the  data  point, 
such  that  Z  =  z  —  zc,  where  z  is  height  above  the  Earth’s 
surface.  As  in  the  study  by  McLandress  et  al  [2000],  we 
specify  the  instantaneous  horizontal  viewing  geometry  in 
terms  of  the  horizontal  azimuth  angle  aT  of  the  orbital 
motion  vector  from  due  east,  and  then  rotate  the  (x,  y,  Z) 
axes  by  this  angle  to  new  axes  (X,  Y,  Z)  such  that  satellite 
motion  now  occurs  along  the  X  axis  and  MLS  views  the 
atmospheric  limb  along  the  Y  axis,  at  90°  to  the  spacecraft 
motion,  as  depicted  in  Figure  6: 


Wz(Z)  =  exp 


(6) 


where  wc  ~  6.8  km  for  channels  3  and  13  [McLandress  et 
al ,  2000],  as  plotted  in  Figure  8b.  The  width  of  the  MLS 
antenna  pattern  yields  additional  smearing  across  the  MLS 
line  of  sight  (LOS)  in  the  Y-Z  plane.  On  ignoring  the  small 
increases  in  this  spatial  width  with  increasing  Y,  this  across- 
LOS  weighting  function  at  the  data  point  can  be 
approximated  as 


WLqs(Y,Z)  =  exp 


z-  m2' 

,  w*  ) 


(7) 


where  by  geometry  the  antenna’s  full  width  half  maximum 
(FWHM)  of  ^0.206°  [Jarnot  et  al ,  1996]  corresponds  to 
Wb  ~  4.85  km  for  channels  3  and  13.  This  function  is 


(X  \  /  cos olt  sin 0 \  /  x  \ 

Y  =  —  sin  olt  cos  0  y  (4) 

Zj  \  0  0  1  )\z-zc) 

[54]  Figure  7  plots  track  angles  aT  as  a  function  of 
latitude  for  both  ascending  and  descending  orbits  and  for 
the  two  yaw  cycles  where  MLS  points  either  to  the  north  or 
to  the  south.  For  south-viewing  geometries,  we  adopt  a 
slightly  different  definition  of  aT  from  McLandress  et  al 
[2000],  such  that  the  values  we  use  are  180°  different.  The 
UARS  yaw  maneuver  is  succinctly  described  by  a  reflection 
of  the  Y axis  about  the  X  axis  (see  Figure  6),  the  convention 
used  by  McLandress  et  al  [2000].  However,  this  has  the 
disadvantage  of  transforming  the  (X,  Y)  axes  from  right- 
handed  to  left-handed  [Anton,  2000]  which  complicates  sign 
conventions  for  GW  wave  numbers.  To  get  around  this,  our 


-90  -60  -30  0  30  60  90 

latitude  (°N) 


Figure  7.  Track  angles  aT  as  a  function  of  latitude  for 
both  ascending  (solid  lines)  and  descending  orbits  (dashed 
lines).  North-  and  south- viewing  orbital  geometries  are 
shown  with  black  and  grey  curves,  respectively.  The  satellite 
motion  along  each  orbit  segment  is  shown  with  arrows. 
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Figure  8.  (a)  MLS  LOS  at  an  angle  (3  below  the  horizontal 
and  the  symmetric  response  around  this  LOS  due  to  the 
antenna  width;  (b)  Gaussian  fit  to  the  vertical  temperature 
weighting  function  for  MLS  channels  3  and  13;  (c)  net  2-D 
MLS  weighting  function  WYz  (f  Z). 


plotted  in  Figure  8a  and  shows  the  MLS  LOS  at  an  angle  (3 
below  the  horizontal  and  the  symmetric  response  around 
this  LOS  due  to  the  antenna  width.  The  net  weighting 
function  along  the  ( Y,  Z)  plane,  centered  about  the  data 
point,  is  the  product  of  equations  (6)  and  (7): 

Wrz(Y,Z)  =  Wz(Z)WLOS(r,Z )  (8) 

which  is  plotted  in  Figure  8c. 

[57]  The  along-track  ( X)  weighting  function  is  determined 
by  a  somewhat  broader  antenna  width  of  ^0.43°  (FWHM) 
in  this  direction  [Jarnot  et  al,  1996],  which  can  be 
approximated  as 


Wx(X)  =  exp 


(9) 


where  «  10.1  km  for  channels  3  and  13  [McLandress  et 
al ,  2000].  Not  considered  in  equation  (9)  is  an  additional 
along-track  broadening  due  to  the  2  s  integration  time  for 
each  MLS  measurement,  which,  given  the  7.5  km  s-1 


satellite  velocity,  corresponds  to  a  15  km  footprint  along- 
track.  We  modeled  this  effect  numerically,  and  found  that  it 
let  to  a  broadened  normalized  Gaussian  that  was  well  fitted 
by  equation  (9)  but  using  «  12  km. 

[58]  Thus  the  three-dimensional  (3-D)  MLS  weighting 
function  about  the  data  point  is 

W(X,YjZ)  =  Wx(X)Wyz(Y,Z )  (10) 

[59]  Equations  (6) -(8)  can  be  re-expressed  by  rotating 
axes  by  an  angle 


2(\>  =  tan  1 


2  [3 

.1  -  [32  +  y~'_ 


such  that 


X' 

r 

z' 


0 

cos<^> 
—  siru\) 


(ii) 


(12) 


This  transformation  enables  equation  (8)  to  be  re-expressed 
as  a  separable  two-dimensional  (2-D)  weighting  function 
[McLandress  et  al ,  2000] 


Wrz(Y,Z) 


Wrz'(Y',Z')  =  exp 


(13) 


where 


wc 

wY'  — - ^ - 

[sin2  4>  +  7(sin  —  (3  cos  4>)] 
Wc 

WZ'  —  - 

[cos2  <\)  +  7 (cos  4>  +  (3  sin  cj>)] 


(14) 


[60]  The  cigar-shaped  MLS  weighting  function  in 
Figure  8c  has  a  major  axis  that  is  aligned  at  an  “effective 
line  of  sight”  angle  4>  that  is  more  horizontal  than  the  actual 
MLS  LOS  angle  (3,  as  was  noted  in  earlier  studies  (e.g., 
Figure  6b  of  Wu  and  Waters  [1997]).  The  angle  4>  quantifies 
this  effective  LOS:  4>  ~  —2.94°  for  channels  3  and  13. 

[61]  The  separable  3-D  MLS  weighting  function  W'(X,  Y', 
Z')  =  Wy(X)  WyZ'  (Y,  Zr)  represents  a  point  spread  function 
about  the  atmospheric  data  point  to  be  measured.  This  point 
spread  function  (or  instrument  function)  will  be  convolved 
through  the  spatially  varying  atmospheric  temperature  field, 
nx,  r,  z\  as  the  satellite  moves  along  the  X  axis  acquiring 
atmospheric  radiances  along  the  limb  track.  The  Fourier 
transform  of  W'  (X,  Y/,  Zr)  is 


R(kx,  ky ,  kz>)  —  Rx (kx)Ryf  {ky  )RZ'  (kz! ) 

=  exp^-('KWxkx)2-('KWyky)2-('KWZ'kZ')2^  (15) 

where  kx,  kY<  and  kz>  are  wave  numbers  (in  cycles  km-1) 
along  the  X,  Y'  and  Z!  axes,  respectively.  Since  convolution 
of  W1  (X,  Y\  Z')  with  L(X,  Y,  Z')  in  the  spatial  domain 
corresponds  the  multiplication  of  the  same  quantities  in  the 
Fourier  domain,  then  R(kx,  kr,  kz)  represents  the  spectral 
sensitivity  or  “visibility”  [Kraus,  1986]  of  the  MLS 


13  of  22 


D03107 


JIANG  ET  AL.:  MOUNTAIN  WAVES  IN  THE  NORTHERN  WINTER  STRATOSPHERE 


D03107 


200  km  500  km  oo  500  km  200  km 


Figure  9.  (a)  MLS  along-track  visibility  Rx  fe)  (solid 
dark  curve)  and  the  band-pass  portion  between  cutoff 
wavelengths  of  30  km  and  90  km  (light  solid  curve).  The 
black  dash-dotted  curve  shows  the  modified  along-track 
visibility  Rrx{k^)  from  equation  (21).  (b)  Two-dimensional 
MLS  spectral  visibility  function  Rr  (kY)  Rz  (kzd-  Grey  lines 
are  explained  in  the  text. 


instrument  to  temperature  structure  of  a  spatial  scale  defined 
by  the  three-dimensional  wave  number  (kx,  kr,  kz),  such  that 

(kx\  ( 10  0  \  (kx\ 

[  kY>  ]  =  0  cos4>  sin  4>  ]  \  kY  ]  (16) 

\  kz>  )  \  o  -  sin  4  cos  <\>  J  \kz  J 

where  kY  and  kz  are  wave  number  components  (in  cycles 
km-1)  along  the  Y  and  Z  axes.  For  a  GW  with  a  wave 
number  vector  K=  (&,  /,  m)  defined  along  the  usual  (x,  y,  z) 
axes,  such  that  Kh  =  (k,l)  =  Kh( coscp,  sincp)  where  Kh  = 
( k 2  +  /2)1/2,  then  kx=  Kh  cos[cp  —  ar],  kY=  Kh  sin[cp  —  ar], 
and  thus 

(kx\  /  cos  a t  —  sin  0  \  /  k  \ 

kY  J  =  I  sin  a T  cos  a T  0  ]  I  /  ]  (17) 

kz)  \  0  0  1  )\m) 

[62]  Thus  equations  (16)  and  (17)  transform  us  from/to 
the  3-D  GW  wave  number  ( k ,  l,  m)  to/from  the  components 
(kx,  kYh  kz)  relevant  to  the  MLS  visibility  calculation  (15). 

[63]  The  sign  of  the  GW  wave  numbers  ( k ,  /,  m)  must 
be  defined  self-consistently  to  properly  represent  via 
equation  (17)  the  precise  3-D  wave  orientations  that 
MLS  is  most  sensitive  to.  The  MWFM  ray-tracing  sign 
convention  for  mountain  waves  described  in  section  4.1 
fulfills  these  requirements. 

[64]  Figure  9b  plots  the  two-dimensional  (2-D)  MLS 
spectral  visibility  function  RYZ'  (kYh  kz)  =  RY  ( kY)Rz  (kz) 
in  (kY,  kz )  space  following  application  of  the  inverse  of  the 
transformation  (16).  The  tilted  nature  follows  directly  from 
the  tilt  of  IV yz  ( Y,  Z)  about  the  Y'  axis,  or  the  tilting  angle  <\> 


in  Figure  8c.  Also  plotted  are  the  lines  kY=  —kz sin<^)  (solid 
grey)  and  kY  =  —kz sin(3  (dotted  grey),  which,  for  a  given 
contour  level,  mark  the  shortest  vertical  wavelength  \z  = 
\kz~l\  and  along-LOS  (or  cross-track)  wavelength  \Y  = 
j kY~ 1 1 ,  respectively,  that  yield  an  MLS  visibility  of  a  given 
level  or  greater  at  some  altitude.  For  example,  on  inspecting 
Figure  9b  we  can  see  that  the  wave  with  the  shortest  cross¬ 
track  horizontal  wavelength  \Y  resolved  at  the  40%  level 
(where  the  dotted  line  intercepts  the  0.4  contour)  is  ^200  km 
with  a  corresponding  vertical  wavelength  \z  of  ^(200  km)/ 

|  sin  (3 1  ^16  km.  The  shortest  vertical  wavelength  resolved  at 
the  40%  level  is  ~  13  km  with  a  corresponding  cross-track 
wavelength  \Fof  ~  1 3 1 sin  <\>\  km  ^250  km. 

[65]  We  see  from  Figure  9b  that  the  peak  response  for 
upward-propagating  waves  ( m  <  0,  hence  kz  <  0)  occurs  for 
waves  with  kY  <  0.  As  noted  by  McLandress  et  al  [2000], 
this  occurs  because  in  this  case  the  GW  phase  lines  are 
quasi-parallel  to  the  MLS  LOS  (3  and  effective  LOS  4>,  and 
thus  the  wave  presents  a  rotated  wave  number  kY  ^  0  to 
MLS,  optimizing  its  visibility  (as  depicted  schematically  for 
view  2  in  Figure  2).  In  other  words,  MLS  is  more  sensitive 
to  GWs  that  are  propagating  in  the  intrinsic  frame  toward 
the  instrument,  i.e.,  in  the  negative  direction  along  the  Y 
axis  depicted  in  Figure  6. 

4.3.  GWs  Resolved  in  MLS  Six-Point 
Limb-Tracking  Variances 

[66]  Given  the  complex  three-dimensional  nature  of  the 
MLS  visibility  function,  we  first  apply  these  analytical 
relations  to  single  monochromatic  plane  waves  to  assess 
in  simple  illustrative  terms  the  types  of  GWs  we  might 
expect  to  resolve  in  the  MLS  six-point  variances.  The 
along-track  spectral  visibility  Rx  (k^)  is  plotted  with  the 
thick  dark  grey  curve  in  Figure  9a.  The  2s  MLS  integration 
time  means  that  data  are  acquired  every  1 5  km  along  track, 
and  so  the  Nyquist  along-track  wavelength  (\x) short  = 
30  km.  The  six-point  data  segments  yield  a  longest  (funda¬ 
mental)  along-track  wavelength  (XxWg  =  90  km.  We 
approximate  waves  resolved  along-track  for  now  by  cutting 
off  the  Rx  (kx)  response  at  (\x)short  and  (\x)iong>  as  shown 
with  the  light  grey  solid  curve  in  Figure  9a. 

[67]  (\x)iong  =  90  km  represents  the  longest  GW 
horizontal  wavelength  that  is  fully  resolvable  in  the  six- 
point  MLS  variances.  A  wave  of  (Xx)/o^  =  90  km  is  only 
visible  to  MLS  when  it  propagates  directly  along  the  X 
axis,  since  propagation  in  any  other  direction  leads  it  to 
project  a  longer  wavelength  component  along  the  X  axis, 
and  thus  it  gets  filtered  out  following  the  long  wave¬ 
length  cutoff  in  Figure  9a. 

[68]  To  study  the  directional  sensitivity  to  waves  in  this 
30-90  km  horizontal  wavelength  band  more  quantitatively, 
we  consider  a  specific  example:  a  GW  of  horizontal  wave¬ 
length  \h  =  Kh1  =  80  km  and  vertical  wavelength  Xz  = 
\m~l\  =  20  km.  We  profile  this  wave’s  MLS  visibility 
(equation  (15))  in  Figure  10  as  a  function  of  latitude  and 
wave  propagation  azimuth  cp  for  north-descending  (ND) 
orbits:  The  track  angles  olt  as  a  function  of  latitude  are 
overplotted  from  Figure  7  with  the  dashed  grey  curve. 

[69]  We  see  that  nonzero  MLS  visibilities  only  occur  for 
propagation  azimuths  cp  near  aT  and  aT  +  1 80°,  indicating 
that  MLS  is  primarily  sensitive  to  this  wave  when  it 
propagates  in  the  along-track  ( X)  direction,  as  then  it 
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Figure  10.  MLS  visibility  for  GWs  with  80  km  horizontal 
and  20  km  vertical  wavelengths  as  a  function  of  latitude  and 
wave  propagation  azimuth  cp  for  north-looking  descending 
(ND)  orbits.  Color  bar  scale  is  linear. 


projects  long  cross-track  wavelengths  \Y  that  are  not 
smeared  out  along  the  limb.  However,  the  peak  response 
(red)  is  offset  from  aT  by  ~12°  in  this  case  for  waves 
traveling  slightly  toward  the  MLS  instrument.  The  reason 
for  this  can  be  gleaned  from  Figure  9b,  which  shows  that 
the  peak  visibility  for  a  given  vertical  wavelength  is  located 
along  the  grey  solid  curve  kY/kz  =  —sin  4>.  Since  kY=  Kh  sin 
(<p  —  aT)  for  a  given  GW,  then  the  peak  response  comes 
from  equating  these  two  equations: 

X/z  sin  4> 

sm(cp  —  aT)  = — - -  (18) 

Xz 


[70]  Figure  10  also  shows  a  sharp  cutoff  in  visibility  from 
values  of  ~  0.2  to  zero  at  angles  cp  —  olt  «  —  30°.  This 
corresponds  to  the  onset  of  the  long  wavelength  along-track 
filter  at  Xx  >  (X x)iong  m  90  km  as  depicted  in  Figure  9a: 
Since  Xx  =  \  cos(cp  —  aT)  and  \h  =  80  km,  Xx  exceeds 
90  km  and  the  wave  is  filtered  out  when  (<p  —  aT)  <  —27°. 

[71]  The  sharp  long  wavelength  cutoff  at  (\x)iong  =  90  km 
in  Figure  9a  is  a  coarse  initial  approximation.  To  delineate 
the  true  nature  of  the  long  wavelength  response,  we  con¬ 
ducted  numerical  simulations  in  which  a  randomly  phased 
GW  of  a  given  along-track  wavelength  was  sampled  at  six 
successive  points  every  15  km,  then  a  six-point  variance 
calculation  exactly  like  that  in  equation  (la)  was  performed. 
In  particular,  the  fitting  and  removal  of  the  linear  trend  in 
equation  (la)  modifies  the  Fourier  content  of  the  data  and 
incompletely  removes  longer  wavelength  GWs.  The  simu¬ 
lated  ratio  of  retrieved  to  actual  variance  from  1000  Monte 
Carlo  simulations  for  each  wavelength  is  plotted  as  a 
function  of  \x  in  Figure  11  (solid  curve).  The  thick  grey 
curve  fits  this  dependence  with  the  function 


/(Xx) 


1  +  (X2 /Xi)^ 
1  +  (Xx/Xi)p 


for  Xz  >  X2 


/(Xx)  =  1  for  Xx  <  X2 


(20) 


where  \\  =  100  km,  X2  =  60  km  and  p  =  4.3.  Thus  we 
modify  the  along-track  component  of  the  3-D  visibility  in 
equation  (15),  Rx  (&x),  to  yield  a  new  along-track  visibility 
function 


R'x^x)  -  Rx{kx)f/\kxl)  (21) 

where  the  square  root  in  equation  (21)  reflects  the  fact  that 
equation  (20)  is  a  variance  response,  whereas  equations  (15) 
and  (21)  refer  to  amplitude  sensitivity.  The  response 
(equation  (21))  is  plotted  as  the  dash-dotted  curve  in 
Figure  9a:  We  note  a  much  more  gradual  roll-off  in  the  long 
wavelength  response  along-track. 


For  channels  3/13  (4>  =  —2.93°),  this  yields  cp  —  aT  «  — 
12°.  Note  from  equation  (18)  that  this  offset  scales  as  X/*/Xz, 
and  thus  the  maximum  absolute  offset  in  the  peak  response 
from  aT  is  defined  by  the  maximum  value  of  X/*/Xz.  Since 
the  visibility  in  Figure  9b  is  <20%  for  all  waves  with  Xz  less 
than  ~  10  km,  we  set  (Xz)min  =10  km  and  note  that  (X/z)max  = 
(\x)iong  cos(cp  —  a^),  where  (Xx)iong  =  90  km.  Substituting 
this  into  equation  (18)  yields 

,  f  \  (Xx)/ong  ^  /1Q\ 

tan(<p  -  ar)max  =  — — -  (19) 

Ozlmin 

Since  (Xx)w/(XZ) mm  =  9,  equation  (19)  yields  a  maximum 
absolute  offset  (cp  —  ar)max  ~  —25°  for  the  peak  wave 
response  in  channels  3/13.  Thus  the  peak  responses  of  all 
the  waves  resolved  in  the  six-point  variances  are  clustered 
at  wave  propagation  azimuths  cp  which  are  less  than  25° 
away  from  the  track  angle  aT.  Note  too  that  for  all  resolved 
waves  the  ratio  X//Xz  <  9  which,  from  the  hydrostatic 
dispersion  relation,  implies  that  waves  resolved  in  MLS 
data  generally  have  high  intrinsic  frequencies  [. Alexander , 
1998;  McLandress  et  al. ,  2000]. 


Figure  11.  Black  curve  shows  mean  normalized  variance 
response  of  a  gravity  wave  of  given  along-track  wavelength 
Xx  extracted  using  the  six-point  variance  calculation 
(equation  (la)).  Each  point  is  the  mean  of  1000  computa¬ 
tions  using  a  randomly  phase  gravity  wave.  Grey  curve  is 
the  analytical  functional  fit/(Xx)  given  by  equation  (20). 
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MIS  "Polar  Diagram"  of  Visibility 
to  Gravity  Waves  of  Given  K|, 


Figure  12.  Schematic  “polar  diagram”  of  MLS  six-point 
limb-tracking  visibility  to  a  GW  of  a  given  wave  number  K 
directed  at  various  azimuths  <p. 


[72]  Figure  12  summarizes  the  results  of  this  section  with 
a  schematic  “polar  diagram”  of  the  response  of  MLS  six- 
point  radiance  variances  to  a  GW  of  a  given  horizontal  wave 
number  directed  at  various  azimuths  <p.  We  note  primary 
sensitivity  to  those  waves  whose  propagation  azimuths  are 
aligned  along-track,  with  an  offset  to  waves  propagating 
toward  the  instrument’s  LOS.  This  sensitivity  is  fixed  to  the 
( X ,  Y)  axes,  which  rotate  as  the  track  angle  aT  varies  with 
latitude  and  yaw  cycle,  as  depicted  in  Figures  7  and  10. 
Going  back  to  Figure  2,  the  optimum  viewing  geometry  lies 
somewhere  between  MLS  view  directions  2  and  3.  Thus,  on 
ND  orbits,  MLS  six-point  variances  (cr^ )  Nf)  resolve  zonally 
propagating  waves  at  the  start  and  end  orbit  latitudes  of 
~80°N  and  ^35°S,  respectively,  peaking  for  waves  with  a 
small  southward  component  to  their  wave  vectors.  At 
midorbit  latitudes  of  ~0°-50°N  where  a T  ^  —60°,  and 
assuming  an  offset  for  the  polar  diagram  peak  of  (<p  —  aT)  « 
—  12°  for  a  \h  =  80  km  and  \z  =  20  km  wave,  the  variances 
are  most  sensitive  to  waves  propagating  at  azimuths  ep  of 
approximately  —70°  and  132°,  the  former  directed  18°  east 
of  due  south,  the  latter  directed  42°  west  of  due  north.  This 
narrow  bidirectional  response  of  MLS  to  GWs  is  implicit  in 
earlier  modeling  results  of  McLandress  et  al  [2000]:  See, 
for  example,  their  Figure  14. 

4.4.  MWFM  Simulations  of  MLS  Variances 

[73]  In  our  MWFM  modeling  of  mountain  waves  in  MLS 
variances,  we  used  temperatures,  horizontal  winds  and 


geopotential  heights  from  UKMO  assimilated  meteorolog¬ 
ical  data  [Swinbank  and  O'Neill ,  1994],  since  they  issued 
uninterrupted  daily  (12Z)  analyses  over  the  duration  of  the 
UARS  mission.  These  fields  are  issued  on  a  global  3.75°  x 
2.5°  longitude-latitude  grid  at  18  standard  pressure  levels 
from  1000  hPa  to  0.4  hPa.  We  ran  MWFM-2  in  “hindcasf  ’ 
mode  through  this  analyzed  atmosphere  over  the  entire 
Northern  Hemisphere  during  every  MLS  observation  day 
for  the  period  1994-1997.  Relevant  mountain  wave  ray 
parameters  returned  at  selected  stratospheric  pressure  level 
are  the  wave  packet’s  geographical  location,  wave  number 
vector  and  its  peak  temperature  amplitude  Ti(p). 

[74]  After  each  hemispheric  MWFM-2  forecast,  we  can 
process  all  the  ray  data  through  an  MLS  GW  visibility 
function.  For  a  given  mountain  wave  ray  number  i  that  has 
propagated  through  the  UKMO  atmosphere  to  a  pressure 
level  p,  where  it  has  a  peak  temperature  amplitude  of  7}(p), 
we  scale  the  amplitude  according  to  the  MLS  GW  visibility 
functions  (15)  and  (21)  to  a  new  MLS-measured  mountain 
wave  temperature  amplitude 

7f“(p)  =  [R'x(kx)Rr>  (kr)Rz'  (**)]  t,(p)  (22) 

Note  that  the  MLS  GW  visibility  term  in  equation  (22)  is  a 
function  of  the  transformed  mountain  wave  wave  number 
(kx,  kyt,  kzl),  which,  from  equations  (5),  (6),  (11),  (14),  (16), 
and  (17),  makes  it  also  an  implicit  function  of  observing 
latitude  (MLS  LOS),  yaw  cycle,  ascending  or  descending 
orbit  and  channel  number.  We  concentrate  primarily  on 
channels  3  and  13  (zc  «  38  km)  hereafter. 

[75]  The  raw  ray  data  from  each  hemispheric  MWFM-2 
forecast  are  voluminous  [see,  e.g.,  Eckermann  and  Preusse , 
1999;  Hertzog  et  al ,  2002].  Since  we  must  generate  and 
store  results  from  many  years  of  such  forecasts,  we  need  to 
reduce  this  data  volume  after  each  forecast  prior  to  storage. 
We  choose  to  average  the  ray  data  from  each  daily  forecast 
to  derive  mean-square  peak  MLS-measured  mountain  wave 
amplitudes  T2{p)  within  2.5°  x  2°  grid  boxes  at  various 
pressure  levels  throughout  the  Northern  Hemisphere.  We 
also  compute  other  potentially  useful  quantities:  for  exam¬ 
ple,  standard  deviations  of  the  RMS  MLS-measured  tem¬ 
perature  amplitude,  for  potential  comparison  with  standard 
deviations  of  MLS  radiances  derived  in  previous  studies 
[McLandress  et  al ,  2000].  This  averaging  procedure  is 
similar  to  the  ones  performed  in  some  previous  MWFM 
studies  [Jiang  et  al ,  2002;  Pierce  et  al ,  2003]. 

5.  Comparison  of  MWFM  Results  With  MLS 
Variances 

5.1.  MWFM  Results  Without  MLS  Filtering 

[76]  Figure  13  plots  maps  of  the  gridbox-averaged 
MWFM  mean-square  peak  mountain  wave  temperature 
amplitude,  T2(p ),  at  5  hPa  (z  ~  37  km),  the  approximate 
altitude  of  MLS  Channel  3/13  data  (zc  ~  38  km).  The  data 
are  further  averaged  in  time  for  simulations  in  all  winter 
months  (December-February)  over  all  MLS  limb-track 
observing  days  during  1994-1997.  However,  no  observa¬ 
tional  filtering  of  any  kind  was  applied  to  these  data:  Thus 
they  represent  the  full  (control)  MWFM  mountain  wave 
temperature  climatology  during  the  MLS  limb-track  obser¬ 
vation  period. 
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Figure  13.  MWFM-simulated  T2(p ),  the  mountain  wave  mean-square  peak  temperature  amplitudes 
(K2)  at  5  hPa  (z  ~  37  km)  without  any  observational  filtering  applied.  Results  are  averaged  in  10°  x  4° 
grid  boxes  for  winter  (December- January  February),  years  1994-1997,  and  are  averaged  for  selected 
MLS  limb-track  days.  The  dotted  contour  lines  are  UKMO  assimilated  total  stratospheric  winds  (ms-1) 
for  this  level  averaged  on  the  same  longitude-latitude  grids  for  the  same  MLS  limb-track  days.  The  color 
bar  scale  is  linear. 


[77]  The  T2(p )  maps  in  Figure  13  show  some  similarities 
with  the  MLS  ND  radiance  variances  (d a)nd  in  Figure  la. 
Obvious  areas  of  some  correspondence  are  variance 
enhancements  over  Scandinavia,  southern  Greenland,  Ice¬ 
land,  and  the  northern  Urals,  immediately  suggesting  these 
MLS  variance  enhancements  may  be  mountain  wave-in¬ 
duced.  However,  in  general  there  are  quite  significant 
differences  between  the  MLS  ND  variances  and  T2{p). 
Some  of  the  main  discrepancies  are  as  follows:  (1)  MWFM 
simulates  enhanced  mountain  wave  activity  all  along  the 
west  coast  of  North  America,  from  southern  Alaska  all  the 
way  down  to  Mexico.  MLS  ND  variances  ( ®a)nd  show 
much  less  activity  here.  (2)  (a2 ) ND  in  Figure  la  is  enhanced 
over  the  Queen  Elizabeth  Islands  and  northern  Alaska, 
where  much  less  MWFM  activity's  evident  in  Figure  13. 
(3)  MLS  (ct^)^  and  MWFM  T2{p)  fields  both  show 
enhancements  across  continental  Europe,  but  the  geograph¬ 
ical  distributions  are  different.  For  example,  the  MWFM 
T2{p)  values  typically  maximize  ~  10-20°  further  south 
than  (d a)nd-  (4)  There  is  significant  MLS^ND  variance 
over  northeastern  Europe,  where  MWFM  f2(p)  is  much 
weaker. 

[78]  One  interesting  feature  of  Figure  13  is  the  almost 
total  absence  of  significant  mean  mountain  wave  tempera¬ 
ture  amplitudes  in  the  northern  equatorial  regions  (0- 
15°N),  despite  the  presence  of  some  significant  mountains 
there  (see  Figure  lb). 

5.2.  MWFM  Results  Using  MLS  ND  Filtering 

[79]  To  model  (b2 )  ND  (reproduced  in  Figure  14a  to  aid 
visual  comparisons),  Figure  14b  plots  5  hPa  MWFM 
T^Dip)  fields:  These  are  the  MWFM-generated  mean- 
square  peak  mountain  wave  temperature  amplitudes  after 


application  of  the  MLS  ND  GW  visibility  filter  (22). 
Substantial  departures  from  the  unfiltered  MWFM  control 
climatology  in  Figure  13  are  evident.  The  mosL  obvious 
change  is  in  overall  scale:  f^D{p)  is  smaller  than  T2{p)  by  a 
factor  of  ~103.  This  highlights  how  most  of  the  simulated 
waves  in  Figure  13  are  not  resolved  by  MLS. 

[80]  The  geographical  distributions  differ  considerably 
too,  with  f2n(p)  showing  enhanced  variance  at  higher 
latitudes  than  f2(p ),  with  little  variance  evident  at  latitudes 
equatorward  of  ^30°N.  The  specific  variance  structure  in 
Figure  14b  agrees  much  better  with  the  (oa)nd  data.  The 
four  specific  inconsistencies  that  were  noted  between 
Figures  13  and  la/ 14a  in  section  5.1  are  all  improved  or 
eliminated:  (1)  The  enhanced  MWFM  mountain  wave 
activity  all  along  the  west  coast  of  North  America  in 
Figure  13  is  substantially  suppressed  after  ND  filtering  in 
Figure  14b,  and  now  agrees  much  better  with  the  (d a)nd 
data  in  Figure  14a.  (2)  Greater  MWFM  activity  occurs  in 
Figure  14b  over  northern  Alaska  and  Baffin  Island,  in  better 
agreement  with  the  {v2A)Nn  data.  (3)  Geographical  variance 
enhancements  in  f2ND(p)  across  continental  Europe  agree 
better  with  Figure  14a,  with  MWFM  variance  enhancements 
generally  displaced  northward.  (4)  Significant  MWFM 
f2n(p)  now  exists  over  northeastern  Europe  in  Figure  14b, 
as  observed  in  ( ®a)nd  *n  Figure  14a. 

[si]  Some  differences  still  exist  between  Figures  14a  and 
14b.  MWFM  activity  over  the  Queen  Elizabeth  Islands  is 
“spottier”  than  the  extended  variance  enhancement  ob¬ 
served  in  Figure  14a.  Modeled  mountain  wave  activity  over 
the  Alps  and  Atlas  Mountains  are  larger  than  observed,  and 
there  are  nearby  MLS  variance  enhancements  in  Figure  14a 
away  from  mountains  that  are  not  seen  in  Figure  14b.  There 
is  also  enhanced  ( ®a)nd  over  Notre  Dame  Mountains 
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Figure  14.  (a)  MLS  ND  limb-track  radiance  variances  <j2a  from  channels  3/13  (z  ~  38  km),  from 
Figure  la.  (b)  T^D{p),  the  gridbox-averaged  MWFM  mean-square  peak  mountain  wave  temperature 
amplitudes  at  5  hPa  (z  ^  37  km)  after  application  of  the  MLS  ND  GW  visibility  filter  (equation  (22)). 
Data  in  both  maps  have  been  averaged  on  a  10°  x  4°  longitude-latitude  grid  during  winter  months 
(December-February)  of  1994-1997.  Dotted  contour  are  UKMO  assimilated  mean  total  stratospheric 
winds  (m  s-1)  at  the  indicated  level  during  limb^frack  days.  Both  color  bar  scales  are  linear:  While  both 
color  bar  scales  have  similar  ranges,  (cr a)nd  and  f2(p)?  though  related,  are  essentially  different  and  so  the 
similarity  of  their  absolute  values  is  largely  serendipitous.  For  example,  different  choices  of  IK  and  1^ 
would  have  changed  the  absolute  values  of  T2(p).  Relative  differences  in  each  case,  however,  are 
comparable  and  hence  significant. 


whereas  MWFM  shows  weaker  activity  that  extends  down 
along  the  Appalachians. 

[82]  Nonetheless,  on  the  hemispheric  scale  there  is  notice¬ 
able  similarity  between  the  MWFM  and  MLS  maps  in 
Figure  14,  with  the  comparison  much  improved  over  the 
unfiltered  control  run  in  Figure  13.  Figure  14  leads  us  to 
conclude  that  there  is  considerable  mountain  wave  informa¬ 
tion  buried  within  these  limb-track  MLS  ND  radiance  var¬ 
iances.  Combined  with  the  analysis  results  in  sections  3.3.1 


and  3.3.2,  we  feel  we  can  confidently  attribute  variance 
enhancements  over  three  of  our  four  focus  regions  to  strato¬ 
spheric  mountain  waves:  Scandinavia,  central  Eurasia,  and 
Greenland.  Enhancements  over  the  Queen  Elizabeth  Islands 
seem  to  be  at  least  partly  mountain  wave-related  on  the  basis 
of  our  analysis  and  modeling,  though  MLS  shows  greater 
enhancements  than  were  modeled  by  the  MWFM. 

[83]  Comparing  model  and  data  in  Figure  14,  we  now 
propose  that  the  following  (®a)nd  enhancements  in 
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Figure  15.  Time  series  of  MWFM- simulated  mean-square  peak  mountain  wave  temperature 
amplitudes,  T^D(p)  at  5  hPa  (z  ~  37  km)  over  southern  Greenland  (blue),  central  Eurasia  (red), 
Scandinavia  (green)  and  the  Queen  Elizabeth  Islands  (black)  for  the  period  of  MLS  limb-track 
measurements. 


Figure  14a  are  also  due  to  stratospheric  mountain  waves 
forced  by  flow  over  the  following  mountainous  features, 
which  we  list  from  Figure  lb  on  progressing  from  west  to 
east:  (1)  the  Alps  (~10°E,  45°N);  (2)  northern  Ural  Moun¬ 
tains  (~60°E,  65°N);  (3)  Putoran  Mountains  (^95°E, 
68°N);  (4)  the  Altai  and  Hangay  Mountains  (^90°- 
100°E,  45°-50°N);  (5)  Sayan  Mountains  and  Yablonovyy 
Ranges  (^95-120°E,  55°N);  (6)  Dzhugdzhur,  Stanovoy, 
Dzhagdy  and  Sikhote  Alin  Ranges  (^120°-140°E,  45°- 
60°N);  (7)  Central  and  Koryak  Ranges  (^160°-170°E, 
^55°-60°N);  (8)  Brooks  Range  and  MacKenzie  Mountains 
(~200-230°E,  65°-70°N);  (9)  Colorado  Rockies 
(~250°E,  40°N);  (10)  Iceland  (~340°E,  65°N);  and  (11) 
Northern  Scotland  (~355°E,  55°N). 

[84]  Figure  15  plots  time  series  of  TjjD{p)  at  5  hPa  over 
our  four  mountainous  focus  regions.  In  agreement  with 
Figures  3b  and  3e,  a  strong  annual  cycle  is  simulated  at 
all  4  regions:  Indeed,  essentially  no  MLS  ND-viewable 
MWFM  activity  occurs  outside  of  the  winter  season.  There 
are  hints  that  some  of  the  site-to-site  and  interannual 
variability  evident  in  the  data  is  also  reproduced,  though 
the  sparse  data  makes  statistical  attribution  here  inconclu¬ 
sive.  In  total,  however,  Figure  15  again  supports  a  mountain 
wave  interpretation  for  the  winter  variance  enhancements 
evident  in  Figure  3. 

5.3.  MWFM  Results  With  MLS  NA  Filtering 

[85]  As  has  been  noted  earlier,  NA  MLS  views  are  less 
suited  to  observing  mountain  waves  during  winter.  To 
assess  this  more  quantitatively,  we  performed  MWFM 
simulations  like  those  in  Figure  14b,  but  applied  an  MLS 
north-ascending  (NA)  filter  (equation  (22))  to  the  modeled 
mountain  rays.  Figure  16b  plots  f2NA(p)  at  5hPa,  the  NA 
analogue  of  f2D{p)  in  Figure  14b. 

[86]  We  see  that  the  vast  majority  of  MWFM-simulated 
mountain  wave  activity  is  substantially  suppressed  or  even 
removed  compared  to  the  ND  map  in  Figure  14b,  consistent 
with  the  MLS  data  (Figures  la,  lc,  and  4)  and  our  concept 
that  NA  views  are  less  favorable  for  observing  mountain 
waves  during  northern  winter.  However,  what  little  activity 
remains  in  Figure  16b  shows  impressive  agreement  with  the 


NA  variances  (cbj )  NA  in  Figure  16a.  In  particular,  we  see 
variance  enhancements  over  southern  and  eastern  Green¬ 
land,  Iceland,  the  northern  Urals,  and  Putoran  Mountains  in 
both  model  and  data.  We  even  see  an  enhancement  over  the 
northern  Verkhoyansk  Range  (130°E,  65°-70°N)  in  north¬ 
western  Russia  that  occurs  in  both  MWFM  and  MLS  NA 
variance  maps,  but  was  not  seen  in  the  corresponding  ND 
variances  in  Figure  14.  The  absences  are  just  as  significant 
too:  For  example,  essentially  no  significant  (a a)na  activity 
occurs  over  North  America,  as  simulated  in  the  corres¬ 
ponding  MWFM  f^A(p)  map. 

[87]  While  the  agreement  between  model  and  data  in 
Figure  16  is  not  perfect  in  every  respect,  it  is  far  better  in 
comparison_  to  the  unfiltered  MWFM  mountain  wave  cli¬ 
matology  T2{p)  in  Figure  13,  which  bears  little  resemblance 
to  th z(oa)NA  data.  Figure  16  leads  us  to  the  surprising 
conclusion  that  many  of  the  high-latitude  NA  variance 
enhancements  in  Figures  lc  and  16a  are  in  all  likelihood 
mountain-wave-related. 

[88]  The  tendency  for  these  NA  mountain  wave  signals  to 
cluster  at  high  latitudes  is  due  to  the  almost  northward  MLS 
views  here  at  the  top  of  the  orbit  (a^  «  0°),  where  both  ND 
and  NA  views  should  have  roughly  similar  sensitivity  to 
waves  propagating  east-west  (see  Figures  6, 7,  and  12).  Jiang 
et  al.  [2003]  noted  a  similar  effect  during  south  yaw  cycles 
when  top-of-the-orbit  observations  viewed  variance 
enhancements  near  the  Rockies  with  near-southward  MLS 
viewing.  In  addition,  the  well-known  variability  of  the  winter 
vortex  flow  [e.g.,  Pawson  and  Naujokat ,  1999;  Waugh  and 
Rong ,  2002]  can  potentially  support  stratospheric  mountain 
waves  with  varying  orientations  with  respect  to  this  varying 
flow,  sometimes  yielding  mountain  waves  with  phase  ori¬ 
entations  more  favorable  to  an  NA  view  than  in  more  typical 
undisturbed  eastward  flow  regimes. 

6.  Discussion  and  Summary 

[89]  The  MLS  instrument  on  UARS  remains  arguably  our 
most  important  source  of  global  data  on  gravity  waves  in 
the  middle  atmosphere,  through  the  GW-related  signals 
contained  in  its  along-track  limb  radiances.  With  a  next- 
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(a)  MLS  Limb-Track  NA  Radiance  Variance  Ch  3/13  (38  km):  DJF  1994-97 


Figure  16.  As  for  Figure  14,  but  for  north-ascending  (NA)  MLS  views. 


generation  MLS  due  for  launch  on  the  Earth  Observing 
Satellite  (EOS-Aura)  platform  in  2004  [Waters  et  al .,  1999], 
data  from  these  two  instruments  offer  the  promise  of  an 
invaluable  long-term  database  on  GWs.  Therefore  it  is  very 
important  that  we  develop  a  thorough  understanding  of  the 
GW  content  of  MLS  radiance  data,  to  provide  a  framework 
and  context  for  analysis  and  interpretation  of  newer  data  to 
come. 

[90]  A  puzzling  issue  in  MLS  GW  studies  has  been  the 
lack  of  success  in  finding  mountain  wave  signals  in  MLS 
data  in  the  Northern  Hemisphere.  This  work  largely  resolves 
the  paradox.  Combining  careful  MLS  data  analysis  with 
global  mountain  wave  “hindcast”  modeling  with  the 
MWFM,  we  have  concluded  that  there  is  a  considerable 
amount  of  mountain  wave  information  buried  deep  within 
MLS  limb-track  radiances  at  middle  and  high  northern 
latitudes,  on  both  ND  and  NA  orbits. 


[91]  The  critical  factor  in  our  work  was  a  detailed 
specification  of  how  a  spectrum  of  mountain  waves  with 
different  3-D  wave  numbers  manifests  as  along-track  radi¬ 
ance  perturbations  as  the  MLS  instrument  acquires  limb 
radiances  along  its  orbital  track.  This  work  built  upon 
important  initial  work  on  this  problem  by  Wu  and  Waters 
[1997]  and  McLandress  et  al  [2000].  Our  so-called  “MLS 
GW  visibility  function”  was  derived,  tested  and  then 
implemented  as  a  general  module  within  the  MWFM  to 
simulate  and  interpret  MLS  mountain  wave  measurements. 
Without  this  conduit  between  model  and  data,  comparisons 
between  MWFM  predictions  (Figure  13)  and  observations 
(Figures  la  and  lc)  were  mediocre  at  best,  with  large  areas 
of  ambiguity.  These  areas  of  discrepancy  essentially  repro¬ 
duce  the  difficulties  that  earlier  studies  encountered  in 
trying  to  associate  structure  in  maps  of  MLS  radiance 
variances  from  the  northern  winter  stratosphere  with  moun- 
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tain  waves  [e.g.,  McLandress  et  al. ,  2000;  Jiang  and  Wu , 
2001].  On  application  of  MLS  GW  visibility  functions, 
much  more  conclusive  comparisons  between  model  and 
data  arose,  for  both  ND  and  NA  orbit  measurements 
(Figures  14  and  16). 

[92]  Indeed,  our  model-data  comparisons  led  us  to  asso¬ 
ciate  a  number  of  variance  enhancements  in  MLS  data  with 
mountain  waves  generated  by  flow  over  variety  of  mountain 
ranges  in  the  Northern  Hemisphere.  Three  of  our  four  focus 
regions  (Scandinavia,  southern  Greenland,  central  Eurasia) 
showed  properties  consistent  with  stratospheric  mountain 
wave  enhancements.  Indeed,  mountain  waves  over  south¬ 
eastern  Greenland  appeared  in  both  ND  and  NA  views  in 
both  data  and  model  variance  maps.  The  fourth  site  (Queen 
Elizabeth  Islands)  also  showed  variance  enhancements 
consistent  with  mountain  waves,  but  modeled  activity  here 
was  less  geographically  extensive  than  observed.  MLS 
variance  maps  showed  enhancements  over  a  number  of 
other  mountain  ranges  that  appeared  to  be  mountain 
wave-induced  also,  as  discussed  and  listed  in  sections  5.2 
and  5.3. 

[93]  The  progress  we  have  made  here  raises  the  question: 
Can  MLS  data  now  be  used  to  refine  and  constrain  global 
mountain  wave  models  and  parameterization  schemes,  like 
the  MWFM?  At  present,  the  answer  to  this  is  “possibly.” 
Since  the  MLS  GW  visibility  function  is  so  critical  to 
attaining  good  model-data  comparisons,  any  uncertainties 
in  it  translate  into  ambiguities  in  model-data  comparisons, 
and  hence  in  model  validation.  The  current  MLS  GW 
visibility  function  is  an  analytical  approximation  specific 
to  the  limb-track  data.  Further  numerical  studies  of  MLS 
GW  filtering  and  extensions  to  limb-scan  data  may  be 
needed,  for  both  MLS  UARS  and  the  new  upcoming 
MLS  on  EOS-Aura.  Because  of  these  uncertainties,  specific 
regions  where  model-data  comparisons  are  poor  cannot  at 
present  be  unambiguously  attributed  to  shortcomings  in  the 
model.  Further  work  may  reduce  these  uncertainties.  For 
example,  detailed  site  studies  (e.g.,  at  our  focus  regions)  that 
locate  and  analyze  large  wave  events  may  be  helpful,  to 
assess  whether  detailed  models  with  MLS  visibility  effects 
included  can  reproduce  the  observations  (see,  e.g.,  the 
CRISTA  case  study  of  Preusse  et  al  [2002]). 

[94]  Currently,  though,  discrepancies  between  model  and 
data  in  Figures  14  and  16  may  reflect  several  things.  First, 
as  just  discussed,  they  may  reflect  shortcomings  in  our 
current  specification  of  MLS  GW  visibility  effects.  Second, 
they  may  reflect  episodic  noise  or  errors  in  MLS  data.  We 
believe,  given  the  extensive  averaging  and  quality  control 
on  MLS  variances,  that  error/noise  contamination  effects  are 
uncommon.  Third,  they  may  highlight  shortcomings  in  the 
model  output,  both  in  the  MWFM’s  basic  physics  and/or  in 
the  analyzed  wind  and  temperature  fields  used  to  perform 
the  hindcasts.  Such  defects  are  likely,  and  MLS  data  offer  an 
attractive  global  data  source  to  improve  and  validate  the 
model. 

[95]  A  fourth  possibility,  though,  is  that  model-data 
differences  here  indicate  variance  enhancements  due  to 
GWs  from  nonorographic  sources,  such  as  convective 
storms,  unbalanced  jet  streams  and  fronts  [e.g.,  Fritts  and 
Nastrom ,  1992].  For  example,  the  NA  variance  enhance¬ 
ment  over  the  Pacific  in  Figure  16a  clearly  cannot  be 
mountain  wave-induced.  Over  land,  ND  variance  enhance¬ 


ments  in  Figure  14a  arise  in  a  zone  ~30°-60°E,  50°-55°N 
over  fairly  flat  regions  of  western  Russia  to  the  southern¬ 
most  Urals  in  the  east.  While  some  of  the  easternmost 
activity  here  may  be  mountain  waves  associated  with  flow 
over  the  southern  Urals,  close  correlations  with  MWFM 
data  in  Figure  14b  are  not  evident.  Substantial  winter  storms 
were  documented  during  the  early  1990s  near  the  southern 
Ural  mountains  [Lassig  and  Moealov ,  2000],  as  well  as 
during  February  1993  and  December  1995  (R.  Lassig, 
personal  communication,  2002),  which  both  fall  within 
the  MLS  limb-scan  and  limb-track  dates.  Severe  storms 
can  be  a  strong  source  of  gravity  waves  [e.g.,  Fritts  and 
Alexander ,  2003]  and  could  potentially  explain  some  of  the 
observed  variance  structure  here.  Further  analysis  is  needed 
to  assess  these  and  other  speculations. 

[96]  The  next  generation  MLS  on  EOS-Aura  significantly 
improves  the  vertical  resolution  over  UARS  MLS.  The 
118  GHz  radiometer  has  a  6.5  km  vertical  and  13  km 
horizontal  field  of  view,  which  should  be  able  to  resolve 
GWs  with  vertical  wavelengths  of  ~4  km  and  cross-track 
horizontal  wavelengths  down  to  ^200  km.  Our  analysis 
here  suggests  that  data  from  this  instrument  should  contrib¬ 
ute  significant  observations  on  mountain  waves  in  the 
extratropical  winter  northern  stratosphere. 
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